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ABSTRACT 


A  theoretical /experimental  study  was  made  of  the  rapid  expansion 
of  partially  dissociated  air  in  a  high  enthalpy  wind  tunnel.  The  theo¬ 
retical  model  included  the  important  nitrogen-oxygen  reaction  kinetics, 
as  well  as  vibrational  exchanges.  Pressure  measurements  were  supple¬ 
mented  by  electron-beam  measurements  of  static  (translational)  tem¬ 
perature,  vibrational  temperature,  and  static  density.  The  primary 
result  was  general  confirmation  oi'  finite -rate  flow  calculations. 
Secondarily,  rotational  temperature  determination  by  electron  beam 
was  found  to  require  a  small  empirical  correction,  vibrational- 
relaxation  rates  from  shock  experiments  were  inapplicable  in  rapid 
flow  expansion,  and  the  static  temperature /density  combination  yielded 
the  optimum  comparison  with  theory.  The  nonequilibrium  data  were 
correlated  with  an  entropy  parameter  for  reservoir  pressures  from  5 
to  20  atm  and  temperatures  from  2300  to  5000°  K.  Verification  of  finite- 
rate  theory  at  high  density  was  inferred  from  data  measured  at  large 
area  ratios  plus  the  simple  nature  of  frozen  flows  downstream  of  the 
nozzle  throat. 
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SECTION  I 
INTRODUCTION 


In  the  past  decade,  large  strides  have  been  made  in  the  theory  of 
gas  dynamic  flow  expansions  in  which  changing  composition  caused  by 
high-speed  chemical  reactions  is  a  distinctive  part  of  the  flow.  Impetus 
for  these  developments  has  come  from  the  design  requirements  of  high 
performance  rocket  nozzles  and  high  enthalpy  wind  tunnel  nozzles.  The 
first  problem  encountered  was  the  evaluation  of  the  basic  thermodynamic 
properties  of  high  temperature  gas  in  equilibrium.  These  properties 
have  now  been  worked  out  in  great  detail  for  high  enthalpy  air  (Ref.  1), 
and  they  form  the  fundamental  basis  for  all  nozzle  performance  calcu¬ 
lations  for  advanced  wind  tunnels.  The  second  major  problem  en¬ 
countered  was  that  of  formulating  a  theory  of  flow  expansion  which 
accounted  for  the  characteristic  finite  rate  of  various  chemical  reac¬ 
tions  in  relation  to  the  rate  of  thermodynamic  expansion  in  nozzle  flows. 
Beginning  with  a  model  of  a  simple  binary  gas  subject  to  sudden  freezing 
of  chemical  reactions,  the  finite-rate  expansion  theories  have  evolved 
to  the  current  theories  which  allow  for  an  arbitrary  number  of  gas  con¬ 
stituents,  an  arbitrary  number  of  interacting  chemical  reactions  and 
vibrational  energy  exchanges,  and  gradual  departures  from  equilibrium 
and  approach  to  frozen  flow  (Refs.  2  through  6). 

Nozzle  flows  producing  vibrational,  chemical,  and  ionic  nonequilib¬ 
rium  have  been  studied,  but  generally,  the  flow  regimes  have  been  such  . 
that  the  translational  and  rotational  degrees  of  freedom  have  been 
assumed  to  remain  equilibrated.  It  is  generally  accepted  that  the  latest 
theoretical  formulations  reproduce  in  every  important  respect  the  actual 
physical  and  chemical  phenomena  prevailing  in  nonequilibrium  nozzle 
flows. 

Experimental  verification  has  lagged  behind  the  development  of  this 
nonequilibrium  theory,  mainly  because  of  the  increased  demands  placed 
on  instrumentation  technology.  It  is  necessary  to  distinguish  between 
verification  of  a  flow  theory  and  flow  calibrations  made  subsequent  to 
that  verification.  For  example,  the  thermodynamic  model  of  flow  in  a 
low  temperature  wind  tunnel  nozzle  is  an  isentropic  expansion  with  a 
constant  specific  heat  ratio.  Verification  that  this  model  correctly 
represents  the  flow  process  requires,  in  addition  to  measurement  of 
the  state  in  the  nozzle  reservoir,  measurement  of  two  independent 
variables  at  some  point  in  the  flowing  gas.  Once  the  adequacy  of  the 
flow  model  is  established,  however,  further  routine  flow  determinations 
require  only  a  single  measurement  in  the  flowing  gas,  and  this  is  found 
to  be  generally  true  regardless  of  the  complexity  of  the  flow  model.  The 
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validation  of  the  flow  model,  on  the  other  hand,  definitely  increases  in 
complexity  as  the  flow  model  increases  in  complexity.  In  principle, 
one  more  independent  measurement  in  the  flow  is  required  for  each 
additional  degree  of  freedom  which  is  activated  at  high  temperature  in 
the  reservoir  and  can  potentially  depart  from  equilibrium  in  the  expan¬ 
sion.  Because  the  number  of  independent  measurements  that  can  be 
made  in  a  flowing  gas  is  limited,  exact  verification  of  a  finite-rate 
expansion  theory  is  feasible  only  for  the  simplest  gases.  The  degrees 
of  freedom  activated  in  high  temperature  air  far  outnumber  the  avail¬ 
able  measurements,  so  that  verification  of  a  flow  model  in  the  exact 
sense  is  currently  unattainable.  This  is  not  to  say  that  partial  verifi¬ 
cation  is  not  possible  or  desirable.  Confidence  in  a  nonequilibrium  air 
model  can  be  maximized,  if  not  made  absolute,  by  simply  making  as 
many  flow  measurements  as  possible,  particularly  if  some  of  these 
measurements  are  especially  sensitive  to  nonequilibrium  effects. 

The  first  attempts  to  check  nonequilibrium  theories  in  air  were 
limited  to  the  conventional  pressure  measurements  in  the  flow,  static 
and  pitot  pressures.  Pitot  pressure  is  relatively  insensitive,  but 
static  pressure  is  moderately  sensitive  to  nonequilibrium  effects,  so 
that  at  least  an  initial  confidence  level  was  established  when  these 
measurements  exhibited  reasonable  agreement  with  contemporary 
theories  (Refs.  7  and  8).  In  recent  years,  a  new  group  of  diagnostic 
measurements  has  become  available  with  the  development  of  the 
electron-beam  probe.  In  its  various  forms,  this  technique  has  added 
static  temperature,  static  density,  and  vibrational  temperature  to  the 
list  of  parameters  which  can  be  measured  in  a  flowing  gas.  Although 
it  has  been  most  extensively  applied  in  unheated  low  density  free  jets 
and  nozzle  expansions,  it  has  recently  been  applied  to  a  number  of  high 
enthalpy  nozzle  expansions  in  which  departures  from  equilibrium  occur 
(Refs.  9  and  10).  It  represents  an  especially  valuable  technique  for 
such  studies  because  the  static  (translational  and  rotational)  tempera¬ 
ture  is  the  most  sensitive  thermodynamic  parameter  to  nonequilibrium 
effects  in  a  flowing  gas.  These  increased  measuring  capabilities  have 
already  indicated  needed  corrections  to  finite-rate  theories  by  confirm¬ 
ing  and  strengthening  earlier  evidence  (Ref.  11)  that  vibrational  relaxa¬ 
tion  rates  obtained  in  shock  experiments  are  not  valid  for  rapidly 
expanding  flows  (Refs.  10,  12,  and  13).  This  fact  has  even  more 
recently  been  substantiated  by  an  infrared  band-reversal  technique 
(Ref.  14). 

A  general  outline  of  the  electron-beam  technique  of  temperature 
measurements  as  originally  formulated  by  Muntz  (Ref.  15)  and  later 
modified  by  others  (Refs.  16  and  17)  is  given  in  Appendix  III.  Several 
different  techniques  have  been  used  to  measure  density  by  an  electron 
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beam  (Refs.  10,  15,  and  18),  but  the  scattering  technique  reported  in 
Refs.  12  and  19  was  used  in  the  durrent  investigation  and  is  briefly 
described  in  Appendix  IV.  The  electron  beam  is  a  very  powerful  diag¬ 
nostic  technique,  but  it  is  inherently  limited  to  very  low  gas  pressures, 
normally  1  torr  and  below.  Consequently,  the  finite-rate  expansions 
of  practical  interest  can  be  studied  by  this  technique  only  after  the  ex¬ 
pansion  has  proceeded  to  very  large  area  ratios  and  high  hypersonic 
Mach  numbers,  far  beyond  the  point  in  the  flow  where  the  departures 
from  equilibrium  occur.  It  will  be  seen  later  that,  in  many  respects, 
this  is  not  really  a  very  serious  limitation.  There  is  also  a  limitation 
expressed  by  the  current  uncertainty  over  the  accuracy  of  the  excitation- 
emission  model  of  Ref.  15.  The  results  of  some  experiments  have 
seemed  to  indicate  that,  for  purposes  of  temperature  determination  at 
least,  an  empirical  correction  to  this  model  must  be  made  (Refs.  16 
and  20).  The  rather  considerable  effort  in  this  field  in  the  past  few 
years  is  summarized  in  Appendix  III.  Fortunately,  in  the  static  tem¬ 
perature  range  used  in  the  present  study,  this  uncertainty  is  only 
10  percent  or  less. 

In  the  present  work,  an  effort  was  made  to  extend  the  experimental 
verification  of  a  finite -rate -expansion  theory  to  a  sufficient  range  in 
reservoir  conditions  that  the  correlation  of  the  nonequilibrium  effects 
with  reservoir  entropy  could  be  explored,  as  done  in  Refs.  21  and  22, 
for  purely  theoretical  solutions.  The  usual  pressure  measurements  in 
the  flow  were  supplemented  by  the  electron-beam  measurements  re¬ 
ferred  to  above,  and  major  attention  was  focused  on  the  stream  static 
temperature,  since  it  displayed  the  greatest  sensitivity  to  nonequilib¬ 
rium.  The  pressure  and  temperature  regime  available  in  the  experi¬ 
mental  facility  to  be  used  was  such  that  it  was  necessary  to  develop  a 
machine  calculation  program  for  a  finite -rate  expansion  allowing  for 
departures  from  equilibrium  upstream  of  a  nozzle  throat  and  in  which 
the  chemical  and  vibrational  reactions  simultaneously  proceeded  at 
finite  rates.  To  this  end,  the  air  model  was  slightly  simplified  by 
neglecting  the  minor  constituents,  argon  (Ar),  carbon  dioxide  (CO2), 
and  all  ionic  species. 


SECTION  II 

EXPERIMENTAL  APPARATUS 


2.1  18-IN.  LOW  DENSITY  WIND  TUNNEL 

The  nonequilibrium  expansions  studied  in  the  investigation  were 
produced  in  a  low  density  hypervelocity  wind  tunnel  supplied  with  air 
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heated  to  high  enthalpy  in  a  250-kw,  d-c  arc  heater  (Figs.  1,  2,  and  3, 
Appendix  I).  A  small  stilling  chamber  is  used  downstream  of  the 
heater,  followed  by  a  simple  conical  nozzle  which  expands  the  flow 
into  an  open-jet  test  section.  The  test  section  is  continuously  pumped 
by  a  six-stage  steam  ejector.  The  reservoir  conditions  produced  for 
this  test  were  total  pressures  from  5  to  20  atm,  total  enthalpies  from 
1150  to  4160  Btu/lbm,  and  stagnation  temperatures  from  2300  to  5000°K. 

The  expansion  nozzle  used  had  a  throat  diameter  of  0.  200  in. ,  an 
exit  diameter  of  18  in.,  a  conical  half-angle  of  10  deg,  and  a  geometric 
area  ratio  of  8100.  With  the  above  reservoir  conditions,  the  static 
pressure  at  the  nozzle  exit  was  in  the  range  from  10  x  10~3  to  31  x  10"^ 
torr,  and  the  static  temperature  at  the  nozzle  exit  was  from  50  to  140°K. 
This  rarefied  flow,  at  mass  flows  from  0.01  to  0.05  lbm/sec,  was 
pumped  through  the  test  section  with  all  six  stages  of  the  steam  ejector 
operating. 


2.2  GENERAL  INSTRUMENTATION 

The  general  location  of  diagnostic  measurements  was  selected  as 
the  exit  plane  of  the  nozzle.  Pitot  pressures  were  measured  1  in.  from 
the  nozzle  centerline  in  a  plane  1  in.  downstream  of  the  nozzle  exit. 

The  pitot -pres sure  measurement  used  directly  for  diagnostic  purposes 
was  made  with  an  uncooled,  flat -faced  tube  having  an  outside  diameter 
of  0.  25  in.  Pitot-pressure  profiles  were  obtained  with  a  flat-faced, 
water-cooled  probe,  0.  5  in.  in  diameter.  Comparison  of  pressures 
measured  by  these  two  probes  indicated  that  low  density  effects  on  pitot 
pressure  were  negligible  for  this  test.  Measured  pitot  pressures 
ranged  from  3  to  8  torr. 

The  static  pressure  was  measured  on  the  nozzle  wall,  1  in.  up¬ 
stream  of  the  nozzle  exit.  The  orifice  diameter  was  0.  25  in.  Condi¬ 
tions  at  this  orifice  were  such  as  to  require  a  low  density  correction  of 
the  type  discussed  in  Refs.  23  through  25.  Capacitance -type  transducers 
were  used  to  measure  both  the  pitot  and  static  pressures. 

Nozzle  reservoir  conditions  were  primarily  determined  by  the 
energy-balance  technique,  although  a  few  checks  of  enthalpy  by  the 
sonic  flow  method  were  made.  The  energy  balance  requires  precise 
measurements  of  heater  input  power,  air  mass  flow,  and  energy  re¬ 
moved  by  cooling  up  to  the  throat.  Input  power  is  determined  by  pre¬ 
cise  voltage  and  current  measurements  at  the  heater  electrodes. 

Cooling  power  is  determined  by  water  volume  flow  measurement  with 
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a  calibrated  turbine-type  flowmeter  and  by  the  temperature  rise  meas¬ 
ured  by  a  ten-element  thermopile.  The  air  mass  flow  is  measured  by 
a  sharp-edged  orifice  with  the  pressure  differential  sensed  by  a  strain- 
gage -type  pressure  transducer.  Inlet  temperature  at  the  orifice  is 
measured  by  a  copper-constantan  thermocouple. 


2.3  ELECTRON  BEAMS 

The  electron  beams  were  projected  across  the  nozzle  exit  flow, 

1  in.  downstream  of  the  nozzle  exit  and  intersecting  the  longitudinal 
axis  of  the  nozzle.  Because  of  conflicting  requirements  of  the  tempera¬ 
ture  and  density  measurements,  two  separate  electron-gun  installations 
were  used.  It  was  not  possible  to  test  with  both  electron  beams  mounted 
in  the  test  cell  simultaneously.  Separate  tunnel  entries  were  used  for 
the  temperature  and  density  measurements  and  the  two  sets  of  data 
were  related  by  using  the  pitot  pressure  as  a  transfer  parameter  for  a 
given  reservoir  condition. 

2.3.1  Temperature  Measurement 

The  temperature  measurements  were  made  with  a  converted 
electron-beam  welding  gun  which  produced  a  2-ma,  20,  000-v  beam. 

This  gun  was  mounted  inside  the  test  section  above  the  free  jet  and  gave 
a  vertical  beam  (Fig.  4).  The  optical  detection  system  consisted  of  a 
21-in.  focal  length  collimating  lens  inside  the  test  section,  a  similar 
condensing  lens  outside  the  test  section,  and  0.  75 -meter  diffraction 
grating  spectrometer  located  outside  the  test  section.  This  system  was 
focused  on  a  small  segment  of  the  vertical  electron  beam  at  the  nozzle 
centerline.  The  radiation  detector  used  on  the  spectrometer  was  a 
Hamamatsu  R106  multiplier  phototube  having  an  S-19  spectral  response. 

2.3.2  Density  Measurement 

Density  measurements  were  made  with  a  television  tube  electron- 
gun  element  which  produced  a  50-ma  beam  at  40,  000  v.  This  beam  was 
injected  horizontally  across  the  nozzle  flow  (Fig.  5).  The  detector  for 
the  scattered  electrons  of  the  density  measurement  was  located  inside 
the  test  section,  shielded  by  a  collimator,  and  placed  above  the  hori¬ 
zontal  beam  in  such  a  position  that  electrons  scattered  out  of  the  beam 
at  an  angle  of  45  deg  from  the  beam  were  detected.  The  detector 
element  itself  was  a  silicon  surface-barrier  detector. 
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SECTION  III 
TEST  PROCEDURE 


3.1  GENERAL 

As  in  most  arc  heater  facilities,  it  is  very  difficult  to  obtain  pre¬ 
cise  repetition  of  reservoir  conditions  from  run  to  run  in  the  18-in.  low 
density  wind  tunnel.  Rather  than  attempting  to  attain  such  repetition, 
the  reservoir  pressure  was  held  as  closely  as  possible  to  a  desired  . 
value,  and  the  stagnation  enthalpy  was  allowed  to  seek  its  own  level, 
determined  by  the  heater  efficiency  and  power  settings  for  each  run. 
The  result  was  a  series  of  runs  at  essentially  constant  ptj  (5,  10,  15, 
and  20  atm)  and  with  ranging  over  rather  wide  limits.  The  basic 
flow  diagnostic  data  measured  at  the  nozzle  exit  were  then  plotted  as 
functions  of  Htl  at  constant  p^.  These  experimental  plots  of  p^,  p^, 
Tj,  T  ,  and  pj  were  faired,  and  the  faired  curves  were  used  to  extract 
values  of  these  parameters  at  even  values  of  Ttl  =  2300,  3000,  4000, 

and  5000°K.  These  are  the  values  tabulated  in  Table  I  (Appendix  II). 

It  has  already  been  noted  that  the  temperature  and  density  data  were 
measured  on  different  tunnel  entries  and  that  pitot  pressure  was  used 
as  a  transfer  parameter  for  these  two  sets  of  data. 


3.2  RESERVOIR  CONDITIONS 

The  reservoir  enthalpy  (Ht^)  was  routinely  determined  by  the  energy- 

balance  method.  The  accuracy  of  this  determination  was  evaluated  two 
ways:  First,  a  theoretical  error  analysis  of  all  the  individual  measure¬ 
ments  entering  into  the  energy-balance  calculation  indicated  an  overall 
accuracy  of  ±3  percent  or  a  total  uncertainty  interval  of  6  percent  for 

(Ht^)  .  Second,  a  number  of  check  calculations  of  enthalpy  by  the 

EB 

sonic  flow  method  were  made,  as  discussed  in  the  next  section. 

A  summary  plot  of  individual  reservoir  conditions  obtained  during 
this  investigation  is  given  in  Fig.  6.  Curves  of  constant  compressibility 
factor  (Z)  are  superimposed  on  these  data  to  indicate  the  degree  of  dis¬ 
sociation  varied  from  essentially  0  percent  below  2500°K  to  a  maximum 
of  20  percent  at  5300°K.  The  reservoir  conditions  for  the  points  referred 
to  in  Section  3.  1  are  given  in  Table  I,  in  terms  of  p^,  T-tp  Htp  E,  Z, 
and  the  specific  heat  ratio  corresponding  to  frozen  conditions  at  the 
reservoir  composition  (Tfrjt^ 
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3.3  SONIC  FLOW  ENTHALPY  CALCULATIONS 

Among  the  factors  entering  into  the  enthalpy  calculation  by  the  sonic 
flow  method,  the  throat  area  was  considered  to  possess  the  greatest  un¬ 
certainty,  especially  during  the  latter  stages  of  a  long  run  when  there 
was  an  unknown  amount  of  deposition  on  the  throat.  This  deposition  is 
the  result  of  progressive  heater  electrode  erosion  during  operation, 
which  leaves  a  fine  deposit  on  the  throat  section.  A  routine  cleaning 
operation  by  brushing  and  air  blasting  is  sufficient  to  remove  most  of 
this  deposit  between  runs.  However,  there  is  often  a  very  thin  persist¬ 
ent  layer  remaining  which  can  only  be  removed  by  ultrasonic  cleaning. 
After  such  a  cleaning,  the  throat  dimension  is  precisely  0.  207  in. ,  as 
determined  by  an  optical  comparator.  Removal  of  the  nozzle  throat 
from  the  tunnel  for  either  ultrasonic  cleaning  or  optical  measurement 
of  d*  is  not  convenient  during  extended  tunnel  operation,  and  it  is  simply 
assumed  that  after  the  routine  cleaning  operation  the  throat  diameter  is 
0.  200  ±  0.  004  in. 

In  Fig.  7  is  given  a  comparison  of  (H^)  with  f°r  the 

initial  point  of  16  runs,  where  d*  was  assumed  equal  to  0.  200  ±  0.004  in., 
or,  in  four  cases,  equal  to  a  precise  optically  measured  value  before 
the  run.  The  tolerance  shown  on  most  of  the  sonic  flow  values  corre¬ 
sponds  to  the  tolerance  in  d*.  There  is  an  evident  statistical  scatter  of 
these  values  about  the  calculated  ±3-percent  scatter  band  of  the  (Ht-.) 

EB 

values.  Note  that  three  of  the  four  points  for  which  d*  was  precisely 
measured  fall  within  this  ±3 -percent  band. 

For  any  points  other  than  the  initial  points  in  a  run,  the  throat 
dimension  was  even  more  uncertain  than  0.  200  ±  0.  004  in.  The  steady, 
but  unknown,  deposition  of  eroded  electrode  material  on  the  throat  pro¬ 
duced  a  systematic  divergence  of  the  two  enthalpy  calculations  for  suc¬ 
cessive  points  in  a  run,  when  a  constant  throat  area  was  used  in  the 
sonic  flow  calculation  {Fig.  8).  Of  course,  there  were  differences  in 
time  between  points  from  run  to  run,  and  probably  differences  in  deposi¬ 
tion  rate  also,  but  an  obvious  upward  trend  in  (H+i)  .  is  shown  in  this 

figure,  which  represents  an  average  of  a  number  of  runs.  The  con¬ 
vergence  of  (Hti)gg  and  ^ti^gonjc  f°r  the  initial  points  in  the  runs, 

however,  is  supporting  evidence  for  the  acceptance  of  the  energy- 
balance  calculation  for  all  points  of  a  run,  since  the  throat  area  does 
not  enter  into  this  calculation. 
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SECTION  IV 
TEST  RESULTS 


In  addition  to  the  reservoir  data.  Table  I  also  gives  tabulated  values 
of  the  test  section  diagnostic  measurements  of  vibrational  temperature, 
rotational  temperature,  density,  static  pressure,  and  pitot  pressure. 
The  pressure  and  density  data  are  nondimensionalized  by  reservoir 
values,  and  effective  area  ratios  for  these  data  are  given  as  obtained 
from  finite-rate  expansion  calculations  described  in  Section  V.  These 
data  are  definitely  experimental  data,  but  they  have  been  interpolated  to 
even  values  of  Tt*.  Each  of  the  test  section  measurements  must  be 

examined  in  the  light  of  certain  qualifications  described  below. 


4.1  VIBRATIONAL  TEMPERATURE 

The  vibrational  temperatures  were  determined  from  nitrogen  {N2) 
molecular  band  intensity  ratios  in  the  electron-beam  emission  (Refs.  10 
and  15  and  Appendix  III).  The  measured  temperatures  were  found  to  be 
uniformly  lower  than  values  predicted  by  finite-rate  expansion  calcula¬ 
tions  using  vibrational  relaxation  rates  obtained  in  shock  tube  studies. 
The  differences  can  be  seen  in  Table  I  as  ranging  from  195  to  1334°K, 
depending  on  the  reservoir  condition.  This  same  result  has  been  found 
in  other  facilities  by  different  techniques,  and  it  is  now  generally  ac¬ 
cepted  that  vibrational  relaxation  times  obtained  in  shock  experiments 
are  not  directly  applicable  in  rapid  expansions.  The  magnifications  re¬ 
quired  on  the  relaxation  rates  to  bring  the  calculated  values  of  vibra¬ 
tional  temperature  into  agreement  with  the  presently  measured  values 
range  from  4.  4  to  28  and  are  tabulated  in  Table  I.  It  was  found  that  in¬ 
creasing  the  vibrational  rates  by  these  amounts  had  a  very  small  effect 
on  the  static  (translational)  temperatures  in  the  finite-rate  expansions. 


4.2  ROTATIONAL  TEMPERATURE 

The  rotational  temperatures  were  determined  by  the  distribution  of 
rotational  line  intensities  within  a  certain  molecular  band  radiation  ob¬ 
served  in  the  electron  beam  (Refs.  9,  15,  and  16  and  Appendix  III). 

Again,  the  radiation  used  is  from  N2  species  so  that,  strictly  speaking, 
the  measured  values  are  N2  molecule  rotational  temperatures  in  the 
flow.  Calculations  have  been  made,  however,  which  support  the  assump¬ 
tion  of  rotational  equilibrium.  Thus,  it  is  a  good  assumption  that  the 
N2  molecule  rotational  temperature  is  equal  to  the  translational  tem¬ 
perature  of  all  species. 
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Because  of  an  uncertainty  as  to  the  basic  theory  of  the  rotational 
temperature  measuring  technique  (Appendix  III),  the  direct  experi¬ 
mental  values  of  Tr  are  supplemented  in  Table  I  by  values  which  have 
been  empirically  corrected  by  the  Robben-Talbot  factor  (Ref.  16).  It 
will  be  seen  in  Section  V  that  the  corrected  values  appear  to  be  more 
consistent  with  the  finite-rate  calculations. 


4.3  DENSITY 

The  static  density  in  the  flow  was  determined  in  a  straightforward 
way  by  counting  electrons  scattered  out  of  a  sensitive  volume  in  a  given 
direction  (Appendix  IV).  This  technique  allows  an  extremely  simple 
and  straightforward  static  calibration  technique.  The  scatter  in  the 
static  calibration  curves  was  found  to  be  approximately  10  percent,  so 
that  the  density  measurement  must  be  considered  to  possess  at  least 
this  degree  of  inaccuracy. 


4.4  STATIC  PRESSURE 

The  static  pressure  was  subject  to  two  appreciable  corrections,  a 
geometric  correction  and  a  low  density  correction.  The  geometric  cor¬ 
rection  was  necessitated  by  the  gradients  in  the  flow  between  the  static 
pressure  orifice  on  the  nozzle  wall  and  the  main  diagnostic  point  on  the 
nozzle  centerline.  The  gradients  transverse  to  the  flow  were  fairly 
large,  as  can  be  seen  in  the  pitot-pressure  profiles  of  Fig.  9,  but  the 
gradients  in  the  streamwise  direction  were  found  to  be  small.  The 
procedure  for  inferring  the  geometric  correction  to  the  static  pressure 
involved  (1)  transfer  of  static  pressure  from  the  wall  orifice  to  the  edge 
of  the  boundary  layer  by  the  standard  assumption  of  negligible  trans¬ 
verse  pressure  gradient  in  the  boundary  layer  and  (2)  relating  of  static 
pressure  at  the  edge  of  the  boundary  layer  to  the  static  pressure  at  the 
nozzle  centerline  by  use  of  the  pitot -pres sure  profiles  and  an  assumed 
is  entropy  in  the  frozen  flow  in  the  test  core. 

The  low  density  correction  to  static  pressure  was  of  the  type  de¬ 
scribed  in  Refs.  23  through  25.  This  low  density  orifice  effect  is 
opposite  in  sense  to  the  geometric  correction;  that  is,  it  increases  the 
inferred  true  static  pressure  over  that  sensed  by  the  orifice.  At  the 
highest  Ttj  and  lowest  p-j- the  low  density  correction  predominates, 

and  the  nett  correction  to  the  measured  static  pressure  is  a  22-percent 
increase.  At  the  lowest  T^  and  highest  pt^,  the  geometric  correction 

predominates,  and  the  net  correction  to  the  measured  static  pressure 
is  a  reduction  of  30  percent  (Table  I). 
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4.5  PITOT  PRESSURE 

The  low  density  effect  on  pitot  pressure  was  calculated  to  be  no 
more  than  from  2  to  3  percent  for  the  range  of  conditions  covered  in  the 
present  test.  This  result  was  confirmed  by  the  agreement  of  p*2  data 
taken  with  the  0.  25 -in.  probe  and  the  0.  50-in.  probe.  Therefore,  no 
low  density  corrections  were  made  to  the  pitot  pressure  in  the  process 
of  analyzing  the  diagnostic  data.  Finite-rate  calculation  of  the  flow  be¬ 
hind  the  probe  shock  wave  was  beyond  the  scope  of  the  present  study, 
and  pitot  pressure  was  simply  assumed  to  be  equal  to  0.  92  p^V-^,  a 
good  approximation  for  hypersonic  frozen  flow. 


SECTION  V 
DISCUSSION 


The  fundamental  purpose  of  this  entire  investigation  was  twofold: 
to  study  the  adequacy  of  the  electron-beam  techniques  as  diagnostic 
tools  and  to  investigate  the  accuracy  of  finite -rate  expansion  calcula¬ 
tions  for  nonequilibrium  flows  in  high  enthalpy  nozzles,  specifically  the 
nozzle  of  the  18-in.  low  density  wind  tunnel. 


5.1  THEORETICAL  FINITE-RATE  NOZZLE  FLOW  CALCULATIONS 

The  development  of  theoretical  models  for  calculation  of  flow  be¬ 
havior  in  nonequilibrium  nozzle  expansions  has  progressed  to  the  point 
where  nearly  every  potential  finite- rate  effect  is  allowed  for  (Ref.  26). 
For  use  in  the  present  study,  a  calculational  procedure  was  developed 
which  is  similar  in  principle  to  that  of  Ref.  17  but  which  is  slightly 
simplified  in  several  respects  to  facilitate  ready  availability  of  solu¬ 
tions.  A  modified  air  model  was  assumed  in  which  the  trace  species 
Ar,  CO2,  neon  (Ne),  and  ionized  species  were  neglected.  Nevertheless, 
all  the  important  features  of  the  energetically  important  N2  and  oxygen 
(O2)  dissociation/recombination  reactions  and  the  nitric  oxide  (NO) 
"shuffle  reactions"  were  included,  as  shown  in  Table  II.  Ionization  was 
neglected  because  the  anticipated  stagnation  temperatures  were  less 
than  6000°K.  The  only  chemical  species  (included  were  N2,  N,  O2,  O, 
and  NO.  In  the  final  programmed  calculation,  the  relative  amounts  of 
N2  and  O2  in  the  low  temperature  composition  can  be  varied  as  desired, 
including  the  extremes  of  pure  N2  or  O2.  For  the  calculations  pre¬ 
sented  herein,  the  composition  was  taken  to  consist  of  0.  79124  moles  of 
N2  and  0.  20876  moles  of  O2,  per  mole  of  air.  The  flow  is  assumed  to 
be  one -dimensional,  and  the  calculation  is  carried  out  by  the  method  of 
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Ref.  27,  using  the  mass,  momentum,  and  energy  conservation  equations, 
the  gas  equation  of  state,  the  caloric  equation  of  state,  the  chemical 
specie  production  equations,  and  the  vibrational  energy  exchange  equa¬ 
tions,  all  in  differential  form.  Internal  energy  of  the  gas  is  determined 
from  the  partition  functions  for  each  chemical  specie,  summed  over  all 
species  present.  The  partition  functions  include  terms  for  translational, 
rotational,  vibrational,  and  electronic  energy,  and  a  heat  of  formation 
term  to  account  for  the  energy  required  to  produce  each  specie  from  the 
standard  state  of  the  gas  at  low  temperature.  The  standard  state  con¬ 
sists  only  of  N2  and  O2  molecules.  When. the  flow  is  essentially  in 
thermodynamic  equilibrium,  as  it  is  assumed  to  be  in  the  reservoir,  the 
chemical  composition  is  determined  from  the  basic  dissociation  equa¬ 
tions  for  the  three  molecular  species.  The  equilibrium  constants  for 
these  reactions  are  determined  from  the  partition  functions  of  the  re¬ 
acting  species.  In  the  nonequilibrium  part  of  the  expansion,  the  total 
production  rate  of  each  specie  is  calculated  from  the  dissociation  reac¬ 
tions  and  the  NO  shuffle  reactions.  These  shuffle  reactions  are  some¬ 
what  faster  than  the  dissociation  reactions,  since  they  require  no  third- 
body  catalytic  agent,  and  thus  tend  to  keep  the  composition  closer  to 
equilibrium.  The  reverse  (recombination)  rates  for  the  reactions  in¬ 
cluded  are  given  in  Table  II.  The  dissociation  rates  are  determined 
from  the  recombination  rates  and  the  equilibrium  constants. 

The  usual  procedure  in  performing  calculations  of  this  type  is  to 
specify  the  nozzle  area  ratio  as  some  analytic  function  of  the  axial 
dimension  (x)  and  to  relate  this  to  the  thermodynamic  solution  by  eval¬ 
uating  the  differential  equations  as  functions  of  the  x-coordinate  also. 
This  works  satisfactorily  when  the  flow  conditions  are  such  as  to  delay 
the  departures  from  equilibrium  until  the  flow  is  in  the  supersonic  part 
of  the  nozzle.  However,  when  nonequilibrium  effects  must  be  included 
in  the  subsonic  flow  and  through  the  throat  region,  the  double-valued 
nature  of  the  area-ratio  function  makes  the  solution  very  difficult  by 
this  method.  In  the  present  case,  the  method  of  Ref.  27  was  used,  in 
which  the  static  temperature  was  taken  as  the  independent  variable,  and 
the  flow  was  related  to  the  nozzle  geometry  by  specifying  the  static  pres¬ 
sure  as  a  function  of  the  x-coordinate.  Although  any  analytic  function 
that  is  monotonically  decreasing  will  do,  a  study  of  several  existing 
nozzle  flow  solutions  indicated  that  the  simple  function 


would  closely  fit  the  calculated  pressure  distributions.  Downstream  of 
the  nozzle  throat  (arbitrarily  at  M  =  1.  10),  the  area  ratio  function  is 
specified.  Currently,  tw?o  successive  cubic  equations  are  used.  The 
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constants  in  the  first  equation  are  selected  so  that  the  flow  properties 
and  their  derivatives  are  continuous  at  the  match  point  between  the 
pressure  function  and  the  area-ratio  function.  The  use  of  static  tem¬ 
perature  as  the  independent  variable  is  a  natural  choice  since  the 
partition  functions,  their  derivatives,  the  chemical  production-rate 
parameters,  and  the  vibrational  relaxation -rate  parameters  are  all 
functions  primarily  of  temperature. 

Particular  care  must  be  exercised  at  the  beginning  of  the  nonequi¬ 
librium  solution.  If  the  flow  is  very  nearly  an  equilibrium  flow, 
arithmetic  inaccuracies  may  result  from  the  requirement  of  finding 
small  differences  between  two  large  numbers.  This  problem  is  dis¬ 
cussed  in  detail  in  Ref.  5.  In  addition,  very  small  step  sizes  are  re¬ 
quired  to  maintain  computational  stability  even  when  sufficient  numerical 
accuracy  is  available.  The  result  is  long  computational  times.  Experi¬ 
ence  has  shown  that  an  equilibrium  solution  may  be  carried  out  for  some 
distance  down  the  nozzle  until  the  chemical  production  rates  indicate  a 
significant  lag  in  the  atomic  recombination  as  compared  to  that  required 
to  maintain  an  equilibrium  composition.  Typically,  the  equilibrium 
calculation  can  be  retained  until  the  actual  recombination  rates  are  ap¬ 
proximately  5  percent  less  than  those  required  by  the  equilibrium  expan¬ 
sion  (Ref.  28).  The  method  used  for  selecting  the  starting  point  of  the 
nonequilibrium  solution  is  essentially  that  of  this  reference,  although 
the  initial  conditions  are  taken  to  be  those  from  the  equilibrium  solu¬ 
tion,  as  proposed  by  Emanuel  and  Vincenti  (Ref.  5). 

The  initial  step  size  for  the  nonequilibrium  solution  is  calculated 
to  correspond  to  the  smallest  characteristic  length  of  the  chemical  re¬ 
action  system  (Appendix  III  of  Ref.  5).  Thereafter  in  the  solution,  the 
maximum  allowable  step  size  is  related  to  the  shortest  characteristic 
length  among  those  from  both  the  chemical  and  vibrational  reactions. 

The  relationship  between  characteristic  length  and  maximum  step  size 
to  retain  arithmetic  stability  depends  on  the  method  of  extrapolation 
used.  According  to  Ref.  5,  a  fourth-order  Runge-Kutta  method  is 
stable  for  step  sizes  up  to  5.7  times  the  minimum  characteristic  length. 
For  a  first-order  (or  linear)  extrapolation  method,  the  allowable  ratio 
is  <  2.  0. 

The  method  of  extrapolating  the  differential  equations  now  in  use  is 
a  compromise  between  the  simplest  assumption  of  constant  gradients 
over  a  small  increment  and  the  complex  calculations  of  the  Runge-Kutta 
method.  The  gradient  of  each  dependent  variable  is  assumed  to  vary 
linearly  with  temperature  about  the  latest  point  in  the  calculations.  The 
rate  of  variation  is  determined  by  comparing  the  last  two  calculated 
values  of  the  gradient.  A  mean  value  of  the  gradient  over  the  succeeding 
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step  is  then  estimated,  and  this  mean  value  is  used  to  calculate  the  next 
value  of  the  variable.  The  dependent  variables  thus  extrapolated  in¬ 
clude  the  static  pressure,  the  vibrational  temperatures,  and  the  chem¬ 
ical  specie  concentrations.  The  density,  enthalpy,  and  velocity  are 
then  calculated,  and  a  new  set  of 'gradients  is  evaluated.  For  each  of 
the  extrapolations,  the  estimated  mean  value  of  the  gradient  is  com¬ 
pared  with  the  last  calculated  value.  Significant  variations  indicate 
either  a  large  step  size  or  large  second  derivatives  of  the  variables, 
suggesting  that  the  linear  approximation  may  not  be  accurate  over  that 
large  an  increment  of  the  flow.  If  the  gradient  varies  across  an  incre¬ 
ment  by  more  than  a  specified  upper  limit,  the  step  size  is  reduced  by 
a  specified  factor,  typically  1.4.  If  all  gradient  variations  are  below 
the  upper  limit  for  ten  consecutive  steps,  they  are  compared  to  a  lower 
limit,  which  is  related  to  the  upper  limit  as  suggested  in  Ref.  5.  When 
all  gradient  variations  are  below  the  lower  limit,  the  step  size  is  in¬ 
creased  by  the  same  specified  factor.  By  this  technique  and  the  imposi¬ 
tion  of  the  stability  limits  on  step  size,  the  calculation  is  allowed  to 
progress  as  rapidly  as  possible  while  maintaining  arithmetic  accuracy. 

The  computational  procedure  described  above  has  been  programmed 
in  FORTRAN  IV  (E)  language  for  calculation  on  an  IBM  360/50  computer. 
Flow  solutions  have  been  obtained  over  the  range  of  reservoir  conditions 
pertinent  to  the  experimental  study  reported  herein  for  the  nozzle  geom¬ 
etry  of  the  18-in.  low  density  wind  tunnel,  described  in  Section  2.  1.  In 
these  cases,  the  nonequilibrium  solutions  were  started  at  temperatures 
within  5  to  65°K  of  the  stilling  chamber  temperature,  corresponding  to 
Mach  numbers  between  0.  13  and  0.  43.  Chemical  "freezing"  occurs 
very  rapidly,  being  essentially  complete  at  the  sonic  point.  Vibrational 
energy  relaxation  continues  further  downstream,  depending  on  the  values 
of  the  rate  parameter  used.  For  these  solutions,  the  average  computer 
run  time  was  approximately  2.  5  min.  For  nozzles  of  a  significantly 
larger  scale  or  for  much  higher  reservoir  pressures,  the  flow  would 
remain  near  equilibrium  much  further  downstream.  In  these  cases, 
the  computer  time  would  be  somewhat  longer,  although  not  unreason¬ 
ably  long. 

Two  additional  features  have  been  developed  and  written  in  the 
FORTRAN  IV  (E)  computer  language,  although  they  have  not  as  yet  been 
utilized  in  actual  calculations.  One  of  these  features  a  set  of  coupling 
factors  which  attempts  to  assess  the  interaction  of  vibrational  nonequi¬ 
librium  on  chemical  reaction  rates  and  chemical  nonequilibrium  on  vi¬ 
brational  relaxation  rates.  The  form  of  the  coupling  factors  was  taken 
from  Ref.  29.  The  second  new  feature  consists  of  the  modified  fourth- 
order  Runge-Kutta  extrapolation  technique  as  described  in  Refs.  26,  30, 
and  31.  If  use  of  the  coupling  factors  or  changes  in  the  flow  properties 
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'  and  nozzle  geometry  result  in  unacceptably  long  computer  run  times, 
the  use  of  this  more  advanced  procedure  may  result  in  significant  gains. 
Several  researchers  have  reported  drastic  reductions  in  computational 
time  as  a  result  of  the  increased  accuracy  which  this  technique  pro¬ 
vides,  thus  permitting  larger  step  sizes. 


5.1.1  Results  of  Finite-Rate  Calculations 


The  net  effect  of  non  equilibrium  on  a  nozzle  flow  may  be  considered 
as  the  combined  effect  of  changes  in  final  chemical  composition  and  con¬ 
sequent  changes  in  the  final  values  of  the  thermodynamic  variables. 
Measurement  of  chemical  composition  of  a  flowing  gas  is  very  difficult. 
Mass  spectrometric  techniques  are  under  development  (Ref.  32)  but  are 
not  yet  considered  sufficiently  developed  to  provide  a  comprehensive  ' 
analysis  of  nonequilibrium  flows.  Thus,  attention  is  directed  to  the 
thermodynamic  variables  which  can  be  measured  with  reasonable  accu¬ 
racy.  The  typical  freezing  of  chemical  reactions  and  vibrational  energy 
exchanges  in  a  finite-rate  nozzle  expansion  leaves  the  flow  in  such  a 
condition  that  atomic  species  and  vibrationally  excited  molecules  are 
unnaturally  present  at  low  temperatures.  The  energy  of  recombination, 
which  would  have  been  released  to  the  flow  had  those  atoms  recombined, 
and  the  vibrational  energy,  which  would  have  been  released  had  the 
vibrational  degrees  of  freedom  remained  equilibrated,  represent  a  with¬ 
holding  of  a  certain  amount  of  thermal  energy  from  the  flow,  in  other 
words,  an  "effective"  heat  removal  from  the  flow.  This  heat  removal 
is  the  primary  effect  of  nonequilibrium  on  the  flow.  It  leads  to  a  de¬ 
crease  of  temperature,  pressure,  and  velocity  and  an  increase  in 
density.  A  secondary  effect  of  nonequilibrium  is  the  imposition  of  non¬ 
unity  values  of  the  gas  compressibility  factor  (Z)  at  low  temperatures. 
This  produces  small  changes  in  the  interrelationships  of  p,  p,  and  T, 
but  the  general  effect  of  heat  removal  is  dominant.  It  is  invariably 
found  that  the  largest  effect  of  nonequilibrium  involves  temperature, 
followed  by  pressure  and  density,  in  that  order. 

It  is  possible  to  express  the  cumulative  nonequilibrium  effects  in 
terms  of  changes  in  the  related  values  of  any  two  of  the  thermodynamic 
parameters.  Since  each  of  these  variables  undergoes  changes  of  vary¬ 
ing  degree,  however,  it  is  often  more  revealing  to  isolate  the  effect  on 
a  single  thermodynamic  variable  by  expressing  it  as  a  function  of  area 
ratio  as  the  independent  variable.  Figure  10  gives  the  variation  of 
static  temperature  with  area  ratio  for  expansions  from  ptl  =  10  atm 

and  T-ti  =  4000°K,  comparing  the  finite-rate  solution  with  various  limit¬ 
ing  cases.  In  fundamental  terms,  the  effect  of  the  kinetic  model  is 
simply  to  change  the  rate  at  which  the  static  temperature  decreases  in 
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the  expansion.  This  temperature  decrease  is  the  result  of  conversion 
of  disordered  energy  (chemical  and  thermal)  to  directed  kinetic  energy 
of  flow.  The  number  of  activated  degrees  of  freedom  which  participate 
in  this  conversion  determines  the  rate  at  which  the  temperature  of  any 
one  of  them  decreases.  Thus,  there  are  the  following  cases  in  order 
of  rate  of  decrease  of  temperature: 

1.  Complete  Equilibrium.  All  of  the  activated  high  temperature 
degrees  of  freedom  participate  in  the  energy  conversion, 
giving  the  slowest  temperature  decrease  with  area  ratio  and 
an  effective  specific  heat  ratio  of  1.  269. 

2.  Finite  Rate.  All  of  the  activated  high  temperature  degrees  of 
freedom  participate  for  only  a  short  initial  segment  of  the  flow, 
after  which  they,  one  by  one,  become  incapable  of  keeping  pace 
with  the  rapidly  changing  thermodynamic  parameters.  The 
result  is  a  greater  rate  of  decrease  in  temperature. 

3.  Completely  Frozen  Expansion.  The  high  temperature  degrees 
of  freedom  are  considered  to  be  so  sluggish  that  they  do  not 
contribute  to  energy  conversion  at  all.  They  remain  activated 
at  the  reservoir  temperature  level,  and  the  kinetic  energy  of 
flow  comes  only  from  the  translational  and  rotational  energy  of 
the  molecules.  The  rate  of  decrease  of  temperature  is  further 
increased.  The  expansion  is  isentropic  with  a  constant 

7=1.  44,  determined  by  the  frozen  composition  in  the  reser¬ 
voir  of  85 -percent  molecules  and  15 -percent  atoms  for  the 
pressure  and  temperature  of  Fig.  10. 

4.  Monatomic  Gas.  This  case  is  not  related  to  the  present  air 
model  but  is  presented  simply  to  show  the  ultimate  limiting 
case  of  expansion  with  the  minimum  possible  number  of  degrees 
of  freedom  (three)  feeding  the  energy  conversion.  This,  of 
course,  gives  the  maximum  rate  of  decrease  of  temperature, 
corresponding  to  7  =  1.  667. 

It  is  noted  that  the  nonequilibrium  effect  on  static  temperature  is 
quite  small  in  the  subsonic  part  of  the  expansion,  but  it  grows  rapidly 
in  the  supersonic  portions.  The  finite -rate  effect  has  decreased  the 
temperature  at  the  throat  only  5  percent  below  the  equilibrium  value, 
whereas  at  A /A*  =  2.  0  the  decrease  is  26  percent,  and  this  grows  to  a 
70-percent  reduction  at  A /A*  =  1000.  At  A /A*  >  100,  the  non  equilibrium 
effect  on  static  temperature  is  nearly  constant.  It  would  be  very  difficult 
to  experimentally  ascertain  finite -rate  effects  in  subsonic  (or  even  sonic) 
flow.  Obviously,  there  is  a  fundamental  advantage  to  be  gained  by  study¬ 
ing  nonequilibrium  effects  at  large  area  ratios,  and  this  is  consistent 
with  the  basic  low  density  limitation  of  the  electron-beam  technique. 
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The  nonequilibrium  effect  on  vibrational  temperature  is  much 
larger  than  on  static  temperature.  The  variation  of  the  N2  molecule 
vibrational  temperature  is  also  shown  in  Fig.  10,  and  it  is  seen  that 
this  temperature  freezes  out  at  a  very  high  level,  3433°K,  whereas  the 
static  temperature  drops  to  very  low  values.  This  curve  results  from 
the  finite-rate  calculation  using  vibrational -rate  data  for  N2  obtained 
in  shock  experiments. 

It  is  interesting  to  observe  the  actual  physical  extent  of  the  regions 
in  the  flow  in  which  the  finite -rate  effects  occur.  These  regions  are 
displayed  in  Fig.  11  for  the  reservoir  condition  of  Fig.  10,  together 
with  part  of  a  full-scale  nozzle  contour.  Nonequilibrium  regions  for 
four  chemical  species  and  three  vibrating  molecules  are  given.  These 
regions  are  arbitrarily  defined  by  a  1 -percent  deviation  from  equilib¬ 
rium  and  an  approach  to  within  1  percent  of  final  frozen  values.  All 
departures  from  equilibrium  occur  upstream  of  the  throat.  It  is 
particularly  interesting  that  the  chemical  finite-rate  region  for  N2  lies 
entirely  upstream  of  the  throat.  Atomic  N  is  not  a  significant  com¬ 
ponent  anywhere  in  this  expansion.  The  total  extent  of  nozzle  flow  in 
which  finite-rate  conditions  exist  is  only  1.9  in.  out  of  a  total  nozzle 
length  of  51  in.  Further,  if  only  the  energetically  important  reactions 
are  considered  (O2  dissociation,  N2  vibration),  the  finite-rate  region 
covers  only  0.  75  in.  around  the  nozzle  throat.  The  important  finite- 
rate  effects  are  completed  by  a  supersonic  area  ratio  of  2.  2,  and  the 
finite-rate  temperature  versus  area-ratio  curve  of  Fig.  10  represents 
a  constant -7  frozen  flow  beyond  this  point.  Flow  velocity  for  equilib¬ 
rium  and  finite-rate  expansion  is  also  given  in  Fig.  11. 


5.2  COMPARISON  OF  THEORY  AND  EXPERIMENT 

The  two  experimental  points  in  Fig.  10  for  vibrational  and  static 
temperature  at  pt^  =  10  atm  and  T^  =  4000°K  show  the  general  result 

found  at  all  reservoir  conditions,  that  the  initial  vibrational  tempera¬ 
ture  calculations  were  much  farther  from  the  measured  data  than  were 
the  static  temperatures. 

5.2.1  Vibrational  Temperature  and  Vibrational  Relaxation  Rate 

The  N2  molecule  vibrational  temperatures  were  determined  by 
analysis  of  the  relative  intensities  of  the  (0,  1)  and  (1,  2)  vibrational 
bands  of  the  N2+  first  negative  system  excited  by  the  electron  beam 
(Appendix  III).  The  vibrational -rate  characteristic  of  the  N2  molecule 
which  was  initially  used  in  the  finite-rate  calculation  was  the  relation 
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rp  =  6.338  x  10“‘s  1 


(2) 


where  Tj  is  in  degrees  Kelvin.  This  is  an  empirical  relation,  with 
form  based  on  transition  probability  theory,  representing  a  large  amount 
of  relaxation  data  obtained  in  shock  wave  experiments,  e.  g. ,  those  re¬ 
ported  in  Ref.  11.  In  all  cases,  the  frozen  vibrational  temperatures  pre¬ 
dicted  by  the  finite- rate  theory  using  this  relation  were  hundreds  of 
degrres  higher  than  values  measured  spectroscopically,  as  shown  in 
Tabi*  I.  This  verified  the  result  reported  in  Ref.  11  and  duplicated  re¬ 
sults  found  in  several  other  laboratories  (Refs.  10,  12,  13,  and  14).  It 
now  setms  to  be  well  established,  by  three  different  experimental  tech¬ 
niques,  that  vibrational  rates  obtained  in  shock  experiments  are  not 
representative  of  the  rates  which  prevail  in  rapidly  expanding  flows. 
Agreement  of  calculations  and  measured  vibrational  temperatures  can  be 
forced  by  simply  increasing  the  vibrational  rates  (reducing  the  relaxation 
times  of  Eq.  (2))  by  the  constant  factors  shown  in  Table  I.  These  rate 
magnifications  were  found  to  vary  from  4.  4  to  28,  depending  on  the 
reservoir  condition. 

The  vibrational -relaxation  time  factor  for  N2  (Tp)  is  given  as  a 
function  of  static  (translational)  temperature  in  Fig.  12.  The  curve 
represented  by  Eq.  (2)  is  shown,  as  well  as  some  of  the  experimental 
data  (from  Ref.  11),  by  means  of  which  the  constants  in  Eq.  (2)  were 
derived.  These  are  the  shock-wave  relaxation  data.  The  reduced 
values  of  Tp  necessary  to  force  agreement  of  experiment  with  theory 
in  the  present  case  are  indicated  by  the  shaded  regions.  These  are 
compared  with  similar  values  obtained  by  a  sodium-line  reversal 
experiment  in  a  pure  N2  expansion  (Ref.  11)  and  with  values  obtained 
in  an  experiment  similar  to  the  present  one  (Ref.  13).  The  present 
results  are  in  good  agreement  with  Ref.  11  at  high  stagnation  conditions 
but  are  generally  higher  than  those  of  Ref.  13  at  all  temperatures.  For 
the  reservoir  conditions  of  this  investigation,  the  calculations  showed 
little  or  no  vibrational  relaxation  occurred  at  static  temperatures  below 
1400°K.  Thus,  the  low  temperature  differences  in  assumed  Tp  between 
Refs.  11  and  13  and  the  present  study  have  no  effect  on  the  resulting 
frozen  vibrational  temperatures. 

It  is  noted  that  the  data  of  Ref.  11  were  for  pure  N2,  whereas  the 
present  data  were  for  air.  In  Ref.  13,  data  from  both  N2  and  air  ex¬ 
pansions  were  subjected  to  correlation  by  a  single  analytic  function. 

It  is  shown  theoretically  in  Refs.  33  through  35  that  N2  should  relax 
more  rapidly  in  air  than  in  the  pure  state  because  of  the  presence  of 
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atomic  species  and  the  exchange  of  vibrational  energy  with  O2  and  NO. 
Theoretical  values  of  Tp  for  pure  N2  from  Refs.  33  through  35  are  given 
in  Fig.  12.  Theoretical  values  of  Tp  for  N2  in  an  equilibrium  air  expan¬ 
sion  lie  between  the  correlation  of  Ref.  13  and  the  present  data  for  the 
higher  reservoir  temperatures,  indicating  general  agreement  with  the 
present  experimental  values.  On  the  other  hand,  the  lack  of  agreement 
between  pure  N2  theory  and  the  data  of  Ref.  11  cannot  be  explained  at 
the  present  time. 


5.2.2  Variables  Other  than  Vibrational  Temperature 

The  effect  of  nonequilibrium  in  an  expanding  flow  will  be  manifested 
by  many  changes  in  the  interrelationships  between  various  thermo¬ 
dynamic  variables.  In  the  present  case,  there  are  four  independent 
measured  variables  available  for  comparison  with  theory,  since  the 
vibrational  temperature  has  already  been  used  to  infer  the  proper  vibra¬ 
tional  rate  to  use  in  the  theory.  These  measurements  may  be  compared 
with  theory  in  many  different  ways,  each  having  its  own  advantages  and 
disadvantages.  Two  such  methods  are  discussed  below. 

5.2.2. 1  Graphical  Comparisons 

One  common  method  Compares  one  of  the  experimental  measure¬ 
ments  with  theoretical  values  obtained  at  effective  area  ratios  deter¬ 
mined  by  the  other  measurements.  Usually,  the  primary  variable 
chosen  for  comparison  has  a  large  sensitivity  to  finite-rate  effects,  e.  g. , 
the  static  temperature.  In  Fig.  13a,  b,  and  c,  the  static  temperature  is 
plotted  versus  area  ratio  for  three  different  reservoir  conditions  selected 
to  cover  the  complete  range  studied.  These  plots  are  simply  enlarged 
segments  of  curves  similar  to  those  of  Fig.  10.  Theoretical  finite-rate 
curves  are  given  for  both  the  normal  and  magnified  vibrational  rates, 
showing  the  relatively  small  effect  of  the  vibrational  rate  on  the  static 
temperature.  The  experimental  static  temperature  (N2  rotational  tem¬ 
perature)  is  plotted  at  "effective"  area  ratios  extracted  from  the  finite- 
rate  solution  at  experimental  values  of  static  pressure,  static  density, 
and  pitot  pressure.  The  static  pressure  measurement  was  subjected  to 
the  corrections  discussed  in  Section  4.  4.  The  temperature  data  points 
represent  the  measurements  as  corrected  by  the  empirical  correction 
of  Ref.  16,  and  the  vertical  bars  extend  upward  to  the  uncorrected 
values . 

It  can  be  seen  from  these  three  plots  that  the  magnitude  of  the  finite- 
rate  effect  increases  as  stagnation  temperature  increases  and  that  there 
is  general  agreement  of  measurements  and  calculations  at  each  reser¬ 
voir  condition.  This  method  of  comparison  of  theory  and  experiment  has 
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the  advantage  of  economy  of  presentation  (all  four  measured  parameters 
entering  into  a  single  plot)  and  also  gives  a  physical  "feel"  for  the 
geometry  of  the  flow.  However,  it  has  the  distinct  disadvantage  that 
the  comparison  is  not  a  strict  experimental/theoretical  comparison, 
since  the  abscissas  of  the  experimental  points  are  obtained  by  reference 
to  specific  finite-rate  calculations.  Of  course,  these  points  will  con¬ 
verge  on  the  theory  in  the  limit  of  perfect  accuracy  of  theory  and  meas¬ 
urement,  but  the  "experimental"  points  in  practical  cases  will  contain 
some  degree  of  theoretical  bias  inherent  in  the  effective  area  ratios. 

It  is  possible  to  compare  measurements  with  theory  in  such  a  way 
that  it  is  strictly  experiment  versus  theory.  In  Figs.  14a,  b,  and  c, 
static  temperature  is  plotted  versus  static  density.  In  this  case,  both 
coordinates  of  the  experimental  points  are  strictly  experimental  num¬ 
bers,  corrected  where  necessary  (Section  IV)  but  uninfluenced  by  any 
theoretical  solution.  This  advantage  is  gained,  however,  at  the  ex¬ 
pense  of  economy  of  presentation,  since  only  two  measured  variables 
at  a  time  can  be  compared  with  theory  on  such  a  two-dimensional  basis. 
With  four  measured  variables,  a  complete  comparison  with  theory 
would  require  six  separate  plots  such  as  those  of  Fig.  14  at  each  reser¬ 
voir  condition. 


5.2.2.2  Numerical  Comparison 


To  retain  the  purity  of  strictly  experimental /theoretical  comparison 
and  yet  avoid  the  complication  of  such  a  large  number  of  plots,  the  com¬ 
plete  experimental/theoretical  comparison  will  be  presented  in  numerical 
rather  than  graphical  form.  For  this  purpose,  we  define 


(experimental  value  of  variable  i) 


theoretical  value  of  variable  i, 
at  excerimental  value  of  variable  i 


'exp 


V'theor  ^iexp) 


(3) 


In  the  formation  of  a  given  j,  there  are  three  sources  of  error: 
experimental  error  in  vj,  experimental  error  in  Vj,  and  error  in  con¬ 
verting  from  vjexp  to  vij.^eor*  It  can  be  shown  that,  for  small  errors. 


fj  -  ftheorj,j  (4) 

where  the  partial  derivative  is  obtained  at  constant  entropy  in  the  ex¬ 
pansion  and  the  e 's  are  percent  errors.  If  we  let 
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Subscript  1  refer  to  static  temperature,  Tj, 

Subscript  2  refer  to  static  pressure,  p^. 

Subscript  3  refer  to  static  density,  pj_,  and 

Subscript  4  refer  to  pitot  pressure,  pt2» 

average  values  of  the  derivatives  among  these  variables  can  be  ex 
pressed  as: 


1.0 
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0.40 
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Vj  dv; 

d  log  Vj 
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d  log  Vj 
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Because  of  the  range  in  values  of  Di,  j,  it  is  clear  that  errors  in  meas¬ 
urement  of  some  variables  will  be  contracted  in  the  formation  of  Ri  j, 

whereas  others  will  be  magnified.  Thus,  there  will  likely  be  preferred 
combinations  of  variables  for  making  the  experimental /theoretical  com¬ 
parison. 


In  general,  the  errors,  € ,  will  contain  both  systematic  and  random 
components.  In  an  average  of  many  Hi,  j.  the  random  components  will 

disappear,  and  j  will  take  the  form  of  Eq.  (4),  in  which  the  errors 
are  only  the  systematic  errors.  In  the  present  case,  R^  j's  were  aver¬ 
aged  over  the  ten  reservoir  conditions  of  Table  I,  and  the  resulting 
(R^  j)’s  are  given  in  matrix  form  as: 


®w)  - 


X 

1.040 

0.975 

1.074 

0.880 

X 

0.810 

1.146 

1.075 

1.187 

X 

1.279 

0.840 

0.950 

0.790 

X 

(6) 


This  matrix  then  represents  the  complete  experimental /theoretical  com¬ 
parison  averaged  over  the  ten  reservoir  conditions.  Before  averaging, 
the  elements  symmetrically  disposed  around  the  diagonal  of  the  Ri,  j 
matrix  are  simply  related  to  each  other.  In  general,  however,  the  ele¬ 
ments  of  Rij  j  are  not  so  related,  and  the  complete  comparison  requires 


all  12  components.  It  can  be  seen  that  the  agreement  of  measurements 
with  theory  ranges  from  within  -21  percent  for  Rpt  p  to  within 

2,  1 
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+28  percent  for  Rp^  p^-  Although  the  individual  R^  j  range  over  these 
rather  wide  limits,  the  overall  average  of  Ri,  j  is  close  to  unity, 

(KTj)  =  1.004  (  7) 

As  anticipated  by  observation  of  Eq.  (4)  and  the  matrix  of  derivatives, 

Eq.  (5),  there  are  preferred  combinations  of  variables.  The  values  of 
(Ri~)  -  1  are  not  simply  proportional  to  values  of  j,  because  the  indi¬ 
vidual  errors  are  not  the  same  for  all  variables.  However,  a  certain 
similarity  between  [(Ri,  j)  -  1  and  j  can  be  observed.  Obviously, 

combinations  of  variables  which  use  static  temperature  as  the  primary 

variable  of  comparison,  (Rt-i  •)*  are  very  desirable,  and  this  could 

■*■*!)  _ 

almost  be  predicted  from  the  form  of  j.  Conversely,  the  (Ri,Tj) 

which  use  Tj  to  define  the  theoretical  value  of  some  primary  variable 
of  comparison  are  quite  poor. 

Whether  a  given  combination  attains  a  high  quality  by  virtue  of 
small  values  of  the  errors,  e,  or  by  favorable  values  of  the  derivatives, 
Di^  j,  it  is  naturally  advantageous  to  focus  attention  on  the  best  of  these 

combinations.  The  j  which  use  static  temperature  as  the  primary 
variable  of  comparison  —  that  is,  the  first  row  of  the  Ri,  j  matrix  —  are 
displayed  in  bar  chart  form  as  (Rti  j  -  1)  in  Fig.  15.  Values  of 
(RTi  j  -  1)  from  other  matrixes  are  also  shown,  namely  for  the  case 

where  normal  vibrational  rates  are  used  in  the  theory  and  also  for  the 
case  where  the  uncorrected  measured  temperature  was  used. 

Certain  general  conclusions  concerning  the  agreement  of  measure¬ 
ments  and  theory  can  be  made  by  observation  of  Fig.  15: 

1.  The  uncorrected  static  temperatures  are  badly  in  error  in 
comparison  with  those  corrected  by  the  empirical  correction  of 
Ref.  16. 

2.  The  magnified  vibrational  rates  give  a  generally  better 
experimental /theoretical  agreement  than  do  the  normal  rates. 
This  is  not  true  for  some  of  the  components  of  R^  j  which  are 

not  shown  in  Fig.  15,  but,  on  the  average,  it  is  true  for  the 
entire  matrix.  It  is  sufficiently  significant  here  that  the  magni¬ 
fied  rates  merely  do  not  disagree  with  experiment,  since  the 
overriding  justification  for  the  magnified  rates  is  not  the  static 
(translational)  temperature  but  the  measured  vibrational  tem¬ 
perature. 
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3.  In  absolute  terms,  the  best  correlation  of  experiment  and 

theory  is  given  by  the  ratio  (Rt^  p^)*  2.  5  percent.  It  is  con¬ 
sidered  significant,  in  this  connection,  that  Ti  and  p^  are  the 
only  experimental  variables  which  share  the  two  character¬ 
istics  of  being  measured  at  precisely  the  same  location  in  the 
flow,  which  eliminates  need  for  geometric  corrections,  and 
being  measured  without  introduction  of  gross  aerodynamic 
disturbances  in  the  flow. 


5.3  GENERAL  CORRELATION  OF  FINITE-RATE  EFFECTS  IN  NOZZLE  EXPANSION 

The  experimental /theoretical  comparisons  discussed  in  Section  5.2.2 
are  all  based  on  absolute  values  of  various  parameters,  no  reference 
being  made  to  the  magnitude  of  effects  produced  by  nonequilibrium  in  the 
flow.  In  this  section,  the  comparison  of  theory  and  experiment  is 
examined  from  the  standpoint  of  the  magnitude  of  the  finite-rate  effects, 
which  will  be  specified  by  parameters  of  the  form, 

vi  . 

ViEQ 

where  v^  may  be  either  an  experimental  measurement  or  a  calculated 
value,  and  vjgQ  is  the  value  of  vj  in  an  equilibrium  expansion. 

It  is  also  of  interest  to  examine  the  possibility  of  a  general  correla¬ 
tion  of  the  finite -rate  effect,  rather  than  simply  isolated  presentation 
of  vi/vigQ  at  separate  reservoir  conditions.  It  was  observed  a  number 

of  years  ago  (Refs.  21  and  22)  that  there  were  a  number  of  nonequilib¬ 
rium  flow  phenomena  in  calculated  expansions  which  correlated  very 
simply  with  the  reservoir  entropy,  Stj/R.  A  further  remarkable  cor¬ 
relation  was  recently  reported  in  Ref.  36,  in  which  the  entropy  param¬ 
eter  was  supplemented  by  a  rate  parameter,  tc.  This  rate  parameter 
was  derived  for  a  nonequilibrium  flow  at  constant  velocity,  which  is 
approximated  by  flow  in  a  supersonic  nozzle  downstream  of  the  throat 
region.  This  parameter  was  not  useful  in  the  present  experiment,  how¬ 
ever,  because  reservoir  conditions  were  such  that  the  nonequilibrium 
began  upstream  of  the  nozzle  throat  where  there  were  large  variations 
in  the  flow  velocity  (Fig.  11).  Instead,  the  L  parameter  of  Ref.  37, 
which  contains  both  reservoir  entropy  and  nozzle  scale,  was  used: 

2  =  +  0.4  £«  <8) 

In  this  relation,  r*  is  expressed  in  inches.  In  Fig.  16a,  the  finite-rate 
effect  on  static  temperature,  expressed  as  T^/T^q  ,  is  plotted  versus  L 
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for  the  range  of  reservoir  conditions  covered  in  Table  I.  The  experi¬ 
mental  points  represent  the  static  temperature  measured  at  each 
reservoir  condition  divided  by  static  temperature  in  an  equilibrium 
expansion  at  points  determined  separately  by  measured  values  of  p^,  pp 
and  Pt2*  Also,  the  Tj/Tj^q  ratios  from  the  finite-rate  theory  were 

all  taken  at  constant  density  in  order  that  the  theoretical  solutions  would 
be  reduced  to  a  single  curve  for  comparison  with  the  measurements. 
Naturally,  the  finite-rate  theory  with  magnified  vibrational  rates  and  the 
best  (corrected)  experimental  values  were  used  as  a  basis  for  this  com¬ 
parison. 

A  clearly  defined  trend  of  the  theoretical  finite-rate  effect  can  be 
seen  in  Fig.  16a,  and  the  experimental  data  agree  with  this  trend  over 
the  complete  range  of  the  L  parameter  covered.  The  finite -rate  effect 
appears  to  be  negligible  for  L  <  26.  0,  but  the  theoretical  reduction  in 
static  temperature  caused  by  the  nonequilibrium  withholding  of  chemical 
and  thermal  energy  from  the  flow  reaches  85  percent  at  E  =  37.  0.  Ex¬ 
perimentally,  this  reduction  is  found  to  be  85.  3  to  87.  2  percent,  depend¬ 
ing  on  the  measurement  used  to  define  the  point  in  the  expansion.  In 
terms  of  absolute  temperature  differences,  the  maximum  theoretical 
finite -rate  effect  on  temperature  is  a  780°K  decrease,  and  the  measured 
decrease  is  792°K,  based  on  measured  density.  As  noted  previously, 
the  Ti j  pj  combination  does  not  uniformly  give  the  best  agreement  of 
theory  and  experiment  at  each  individual  reservoir  condition  but  aver¬ 
aged  over  all  reservoir  conditions  it  does. 

In  Fig.  16b,  the  finite-rate  effect  on  static  pressure  -is  presented 
as  a  function  of  E  in  the  same  way  as  for  temperature  in  Fig.  16a. 
Although  there  are  reservations  about  the  accuracy  of  the  (corrected) 
static  pressure  measurements,  it  is  considered  instructive  to  make 
this  correlation  because  the  static  pressure  is  nearly  as  sensitive  to 
finite-rate  flow  effects  as  is  the  static  temperature.  As  with  tempera¬ 
ture,  there  is  a  regular  downward  trend  of  theoretical  Pi/Pi^q  as 

L  increases.  There  is  also  observable  a  fair  agreement  of  the  experi¬ 
mental  measurements  with  this  trend,  at  least  in  the  range  29  <  L  <  34. 
The  relatively  larger  scatter  of  the  Pi/Pi^q  P°ints  defined  by  the  other 

measurements  is,  of  course,  related  to  the  inherent  error -magnifying 
effect  of  the  large  derivative  of  static  pressure  with  respect  to  all  other 
variables,  i.  e.,  the  second  row  of  j,  Eq.  (5). 
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5.4  RESULTS  OF  THEORETICAL/EXPERIMENTAL  AGREEMENT 
5.4.1  General  Inferences 

It  seems  to  have  been  demonstrated  reasonably  well  (Figs.  13 
through  16)  that  the  static  temperature  in  a  nonequilibrium  flow  is  accu¬ 
rately  measurable  by  the  electron-beam  technique  and  that  the  measure¬ 
ments  agree  with  finite-rate  expansion  calculations  to  a  very  satisfactory 
degree.  At  various  stages  of  development  of  finite-rate  expansion 
theories,  experimental  verification  has  been  desired  to  assure  that  the 
chemical-kinetic  formulation  used  was  adequate  under  conditions  of 
rapid  expansions  and  also  that  the  basic  characteristic  rate  constants 
of  various  reactions  were  reasonably  accurate.  The  present  results 
represent  a  further  contribution  to  the  growing  confirmation  of  these 
theories  by  virtue  of,  first,  measurement  of  that  thermodynamic  vari¬ 
able  which  is  most  sensitive  to  nonequilibrium  and,  second,  a  large 
enough  range  in  reservoir  conditions  that  a  general  correlation  of  the 
nonequilibrium  effects  is  demonstrated.  For  quite  some  time,  there 
has  been  no  doubt  that  the  reaction  kinetics  itself  was  properly  repre¬ 
sented,  but  the  accuracy  of  the  basic  rate  data  has  been  under  question. 
The  present  experiments  confirm  the  current  belief  that  vibrational  de- 
excitation  rates  which  prevail  in  rapid  expansions  must  be  determined 
in  such  flows  and  cannot  be  inferred  from  vibrational  excitation  rates 
observed  in  shock-type  flows.  This  fundamental  difference  between 
vibrational  rates  in  compression  and  expansion  flows  was  determined 
specifically  only  for  N2,  but  similar  results  have  been  found  for  CO 
(Ref.  14),  and  it  is  very  likely  that  the  same  situation  prevails  for  O2 
and  NO.  On  the  other  hand,  the  chemical  reaction  rate  data,  which 
have  been  mainly  obtained  in  shock  experiments,  appear  to  be  entirely 
valid  in  rapid  expansions. 

Although  the  present  work  represents  a  positive  step  in  the  verifi¬ 
cation  of  finite-rate  flow  expansion  theory,  it  is  far  from  complete. 
Ideally,  at  least  three  more  features  would  be  required.  First,  verifi¬ 
cation  should  be  extended  to  such  a  range  in  reservoir  conditions  that 
a  complete  transition  from  equilibrium  to  fully  frozen  flow  could  be 
measured.  Second,  the  verification  should  be  completed  in  the  exact 
sense  by  performing  measurements  of  composition  in  the  flowing  gas. 
Third,  a  direct  verification  of  theory  would  be  desirable  in  the  pre¬ 
cise  regions  of  flow  in  which  the  nonequilibrium  disturbance  is  gen¬ 
erated,  i.  e. ,  at  small  area  ratios  near  the  nozzle  throat.  The  first 
two  of  these  ideal  requirements  appear  to  be  quite  possible.  However, 
it  is  doubtful  that  the  third  requirement  will  ever  be  met  because  of  the 
very  large  gradients  of  all  flow  properties  and  the  severe  environment 
near  the  nozzle  throat. 
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The  fact  that  thermodynamic  variables  measurable  at  large  area 
ratios  are  evidently  predictable  by  finite -rate  theory  is  quite  strongly 
suggestive  that  this  theory  is  actually  closely  representative  of  the  flow 
in  those  regions  where  finite-rate  conditions  occur  near  the  throat. 

The  reason  for  this  is  that  the  flow  downstream  of  the  (limited)  finite- 
rate  region  is  at  constant  composition  with  a.  constant  specific  heat 
ratio  defined  by  the  frozen  composition  {7  >  1.  4  because  of  the  frozen 
monatomic  species).  The  thermodynamics  of  such  flows  is  especially 
simple  and  well  known;  the  7=1.4  isentropic  expansion  in  a  conven¬ 
tional  wind  tunnel  nozzle  is  the  most  familiar  example  of  such  a  flow. 
Except  for  a  small  amount  of  curvature  at  A /A*  <  10  (Fig.  10),  the 
curve  of  T  versus  A /A*  is  a  straight  line  on  logarithmic  coordinates 
and  can  be  shown  to  have  a  slope  of  ( 1  -  7).  Extension  of  such  curves 
in  either  direction  over  large  ranges  in  area  ratio  is  very  simple  and 
has  a  sound  theoretical  basis.  Hence,  the  good  agreement  of  measure¬ 
ment  and  theory  at  large  area  ratios  is  inferred  to  indicate  the  flow 
solutions  are  valid  at  the  beginning  of  the  frozen  flow  region  at  low 
area  ratios. 

5.4.2  Routine  Flow  Calibration  of  Nonequilibrium  Flows 

The  direct  proof  that  thermodynamic  variables  at  large  area  ratios 
are  predictable  by  finite-rate  expansion  theory  {and  indirect  inference 
that  the  same  is  true  at  low  area  ratios)  provides  a  basis  for  perform¬ 
ing  routine  calibration  of  high  enthalpy  nozzle  expansions  which  depart 
from  equilibrium.  The  diagnostic  procedure  can  be  entirely  similar 
to  that  used  in  conventional  wind  tunnel  practice  and  will  differ  only  in 
the  theoretical  formulation  of  the  thermodynamics  of  the  expansion  pro¬ 
cess.  The  following  steps  will  be  included: 

1.  Precise  determination  of  the  state  of  the  gas  in  the  stagnation 
reservoir.  If  the  reservoir  is  adequately  sized,  the  gas  will 
be  in  equilibrium  at  this  point,  and  its  state  can  be  determined 
by  measurement  of  any  two  thermodynamic  variables.  Reser¬ 
voir  pressure  is  almost  alwrays  one  of  these  measurements, 
and  the  reservoir  temperature  or  enthalpy  is  the  other.  In 
most  arc-heated  facilities,  this  means  accurate  energy -balance 
measurements . 

2.  Analytical  formulation  of  the  finite-rate  thermodynamic  pro¬ 
cesses  in  the  nozzle  expansion  by  the  theory  described  in 
Section  5.  1.  The  complexity  of  the  thermodynamic  properties 
?t  high  temperature  and  the  addition  of  the  expansion  rate  as 
an  important  flow  parameter  dictate  that  a  separate  calculation 
is  required  for  each  reservoir  condition  and  nozzle  size. 
Nevertheless,  these  calculations  are  the  exact  analogue  of  the 
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7=1.4  compressible  flow  tables  used  in  conventional  wind 
tunnel  work  independent  of  the  nozzle  scale  or  supply  condition. 

3.  Measurement  of  a  single  convenient  parameter  at  some  point 
of  interest  in  the  expanded  flow  to  locate  this  point  in  the 
theoretical  calculation.  The  values  of  the  other  flow  variables 
at  that  point  are  thus  defined.  The  remarks  made  in  Section 
5.  2.  2  concerning  the  desirability  of  various  parameters  for 
determination  of  effective  area  ratio  apply  only  to  the  very 
special  conditions  of  this  investigation,  namely,  a  flow -model 
validation  study  performed  at  very  low  densities.  It  is  not 
foreseen,  for  example,  that  the  static  density  measurement 
will  become  a  routine  measurement  for  flow  calibration  nor  do 
the  nonequilibrium  limitations  on  pitot  pressure  apply  generally 
at  higher  densities.  It  is  probable  that  the  pitot  pressure  and 
the  static  pressure,  will  continue  to  be  the  variables  most  suited 
to  routine  calibration  of  high  enthalpy  flows. 

The  requirement  of  a  separate  and  distinct  finite-rate  flow  calcula¬ 
tion  for  each  specific  nozzle  expansion  is  likely  to  prove  very  burden¬ 
some  in  wind  tunnel  work.  The  work  of  Ref.  36  has  pointed  the  way  to 
a  means  of  systematizing  the  results  of  many  such  calculations  so  that 
something  approaching  the  simplicity  of  constant-y  isentropic  flow  tables 
may  be  possible  even  for  highly  nonequilibrium  flows. 


SECTION  VI 
CONCLUSIONS 

1.  A  systematic  study  of  experimental  measurements  in  high 
enthalpy  nonequilibrium  nozzle  flows  of  air  has  revealed  very 
good  agreement  with  conditions  predicted  by  finite -rate  ex¬ 
pansion  calculations  for  the  specific  nozzle  geometry  and 
reservoir  conditions  in  the  range  from  5  to  20  atm  of  pres¬ 
sure  and  temperatures  from  2300  to  5000°K.  The  theoretical/ 
experimental  comparison  was  primarily  in  terms  of  the  static 
temperature  measured  by  the  electron-beam  probe  technique. 
The  flow  expansion  used  a  slightly  simplified  air  model  and 
allowed  for  finite-rate  chemical  reactions  and  simultaneous 
(but  uncoupled)  finite-rate  vibrational  energy  exchanges. 

2.  Vibrational  temperatures  measured  by  the  electron-beam  probe 
were  hundreds  of  degrees  lower  than  values  predicted  using 
vibrational  rates  obtained  in  shock -tube  experiments.  The 
measured  values  corresponded  to  vibrational  rates  4  to  28  times 
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more  rapid  than  the  shock  tube  rates,  which  is  in  general 
agreement  with  several  other  reported  experiments.  Appar¬ 
ently,  vibrational  deexcitation  rates  to  be  used  in  expansion 
flow  calculations  must  be  obtained  in  experiments  on  expand¬ 
ing  flows. 

3.  Vibrational -rate  magnifications  of  up  to  28x  were  found  to 
have  only  small  effects  on  the  calculated  static  translational 
temperatures  in  a  flow,  but  on  the  average,  these  small  effects 
were  in  the  direction  of  better  agreement  with  the  measure¬ 
ments.  The  generally  good  agreement  of  the  static  tempera¬ 
tures  implies  that  the  chemical  reaction  rate  data  used  in  the 
calculations  were  accurate. 

4.  A  general  correlation  of  the  nonequilibrium  effects  on  static 
temperature  and  pressure  in  the  expansions  was  obtained  for 
the  entire  range  of  operating  conditions,  using  a  correlating 
parameter  which  was  dependent  mainly  on  reservoir  entropy. 

5.  The  excellent  experimental/theoretical  agreement  found  at 
large  area  ratios  was  interpreted  as  indirect  verification  of 
nonequilibrium  flow  solutions  at  small  area  ratios  because  of 
the  inherent  simplicity  of  the  frozen,  constant-7  flow  down¬ 
stream  of  the  regions  in  the  nozzle  in  which  finite-rate  condi¬ 
tions  actually  prevail.  The  finite -rate  theory  can  thus  provide 
the  basis  for  routine  calibration  of  high  enthalpy,  nonequilib¬ 
rium  nozzle  expansions  at  all  area  ratios. 

6.  The  electron-beam  techniques  for  measuring  static  temperature 
and  static  density  were  found  to  be  very  valuable  tools  for  pur¬ 
poses  of  validation  of  the  flow  theory.  In  fact,  for  the  special 
conditions  under  which  this  investigation  was  made,  the  density 
and  temperature  were  found  to  be  the  best  combination  of  ex¬ 
perimental  measurements  for  demonstrating  the  validity  of  the 
theory.  For  purposes  of  routine  flow  calibration,  however,  it 
is  likely  that  the  static  and  pitot  pressures  will  remain  as  the 
most  convenient  measurements. 

7.  Evidence  was  found  in  the  experimental  measurements  that  the 
formal  Muntz  theory  of  electron-beam  excitation  and  emission 
is  not  completely  accurate  and  requires  a  small  empirical  cor¬ 
rection.  In  the  temperature  range  of  interest  in  the  present 
study,  this  correction  is  approximately  8  percent  of  the  meas¬ 
ured  rotational  temperature. 

8.  The  present  experiment  represents  the  first  comprehensive 
test  of  the  density  measurement  by  counting  of  electrons 
scattered  out  of  the  electron  beam.  The  repeatability  obtained 
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was  ±10  percent.  This  technique  is  considered  very  promising 
because  it  does  appear  capable  of  further  refinement,  and  it 
has  the  unique  capability  of  measuring  the  mass  density  inde¬ 
pendent  of  the  degree  of  dissociation  of  the  scattering  particles. 
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ELECTRON  GUN 


Electran-Beam  Installation 
Fig.  4  Electron- 


’  Temperature  Measurement 


Installation 
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b.  Photograph  of  Flow  with  Electron  Beam  Used  for  Temperature  Measurement 

Fig.  4  Concluded 
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DATA  POINT  WITHIN  A  RUN 


Fig.  8  Effect  of  Throat  Contamination  on  Indicated  Sonic  Flow  Enthalpy, 
Assumed  d*  =  0.200  ±  0.004  in. 
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VELOCITY,  tt/s#  c 


Fig.  11  Finite-Rate  Flow  Regions  and  Velocity  in  Nozzle,  pt  =  10  atm,  T 
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TEMPERATURE,  »K 


Fig.  12 


Nitrogen  Vibrational  Relaxation  Time  Factor  versus  Temperature 
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a.  p.  =  20  atm,  T,  =  2300°K 

Fig.  13  Experimental  and  Theoretical  Static  Temperature  versus  Area  Ratio 
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EQUIVALENT  INVISCID  AREA  RATIO 


b.  pti  =  10  otm,  Ttf  =  4Q00°K 
Fig.  13  Continued 
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c.  pt  =  5  atm,  Tf  =  5000  °K 
M  M 

Fig.  13  Concluded 
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STATIC  DENSITY  RATIO,  p)/pii 


o.  p,  =  20  atm,  T,  =  2300°K 
'l  T1 

Fig.  14  Experimental  and  Theoretical  Static  Temperature  versus  Static  Density 
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b.  p,  =  10  atm,  T,  =  4000°K 
1  1 
Fig.  14  Continued 
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Fig.  15  Average  Deviation  of  Measured  Static  Temperature  from  Finite>Rate  Theory  for  Ten 

Reservoir  Conditions,  p.  =  5  to  20  atm,  T,  =  2300  to:5000°K 
T1  T1 
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2 

a.  T  ,/T  i  _  _  versus  2 
■  1  E  Q 

Fig.  16  General  Correlation  of  Finite-Rate  Effect  on  Static  Temperature  and  Static  Pressure 
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O  AT  MEASURED  f, 

A  AT  MEASUREDT, 

□  AT  MEASURED  p, 

- FINITE-RATE  THEORY 


TABLE  I 

SUMMARY  OF  RESERVOIR  CONDITIONS  AND  MEASUREMENTS  AT  NOZZLE  EXIT 


m 

-j 


Reservoir  Conditions 

Vibrational  Temperature 

Static  (Rotational) 
Temperature 

Point 

Ptp 

atm 

Ttp 

•K 

Ht,* 

Btu/lbm 

X 

»fr,  t 

Z 

Rxperimental, 

“K 

Theoretical 
(Normal  Rates) 
“K 

Magnification 

Factor 

Required 

Measured, 

”K 

Corrected, 

°K 

1 

IS 

2300 

1154 

1.0003 

1.403 

29.  47 

1937 

2132 

4.  4x 

58.  3 

51.  4 

2 

20 

2300 

1154 

1.  0002 

1.403 

29.  18 

1777 

2099 

9.  3x 

55.8 

48.0 

3 

10 

3000 

1635 

1.  008 

1.407 

31.  34 

2228 

2701 

11.  Ox 

83.2 

75.  3 

4 

IS 

3000 

1625 

1.  006 

1.406 

30.  79 

2136 

2630 

1  1.  3x 

81.8 

73.9 

5 

20 

3000 

1619 

1.  005 

1.4055 

30.  48 

1983 

2575 

17.  Ox 

76.  1 

68.  4 

6 

5 

4000 

2911 

i.  one 

1.449 

34.  87 

2682 

3587 

22.  5x 

123.3 

114.5 

7 

10 

4000 

2755 

1.075 

1.440 

33.  80 

2592 

3433 

15.  3x 

115.2 

106.  5 

8 

15 

4000 

2676 

1.064 

1.434 

33.21 

2508 

3331 

14.  Ox 

112.7 

104.0 

9 

20 

4000 

2625 

1.  057 

1.431 

32.  80 

2373 

3257 

16.  Ox 

110.8 

102.  1 

10 

5 

5000 

4157 

1.  182 

1.484 

37.  14 

2933 

4267 

28.  Ox 

139.2 

130.  3 

Static  Density 

Static  Pressure 

Pitot  Pressure 

Point 

(pt^ti)uxp 

Inferred 

A/A* 

(pl/ptl)exp 

(Pi/Pti> 

1  1  corr 

Inferred 

A/A* 

<Pl2^ptl)exp 

Inferred 

A/A* 

5.  75  x  10'5 

■ 

1.  73  x  lO-6 

1.  23  x  10"° 

4410 

3. 31  x  10"4 

5000 

5.  30  x  10'5 

1.61  x  10‘6 

1.  12  x  10‘° 

4760 

3.  02  x  10"4 

5510 

7.  23  x  10~b 

2.  1  x  10"6 

1.55  x  10'fi 

3810 

4. 0  x  10"4 

4180 

6.  6  x  10'  •’ 

3710 

1. 96  x  10~6 

1.43  x  10'6 

4000 

3.81  x  10-4 

4340 

6.  2  x  lO'0 

3910 

1.  79  X  10-G 

1. 29  x  10'6 

4410 

3.36  x  10"4 

5040 

17.5  x  10'5 

1410 

2. 73  x  10'6 

2.  85  x  10‘6 

2210 

8.  03  x  10"4 

2040 

B 

8.  3  x  10'b 

2870 

2. 34  x  10"6 

1.  88  x  10_li 

3150 

5.02  x  10"4 

3330 

8 

7,  5  x  10'5 

3150 

2.  13  x  10'6 

1. 62  x  10‘6 

3610 

4.  23  x  10'4 

3980 

9 

7,42  x  10"° 

3120 

2.03  x  10"6 

1.55  x  10'° 

3810 

3.85  x  10"4 

4410 

10 

20. 7  x  10'5 

1280 

2.  85  x  10"e 

3.  47  x  10'° 

1800 

9. 36  x  10"4 

1905 

AEDC-TR-69-66 


AE  DC-T  R-69-66 


TABLE  II 

REACTIONS  INCLUDED  IN  FINITE-RATE  FLOW  MODEL 


Reaction 

Cl 

C2 

c3 

n 

O2  4  O2 

-v- 

20  4  O2 

8.  0  x  1013 

-3/2 

0 

2 

O2  +  O 

20  4  O 

2. 3  x  1014 

-3/2 

0 

2 

O2  +  N2 

— 

20  4  N2 

6. 2  x  109 

-1/2 

0 

2 

O2  +  N 

20  4  N 

3.  0  x  109 

-1/2 

0 

2 

O2  4  NO 

T” 

20  4  NO 

3.  0  x  109 

-1/2 

0 

2 

N2  4  O2 

T” 

2N  4  O2 

1.  1  x  1010 

-1/2 

0 

2 

N2  4  O 

2N  4  O 

1.  1  x  10l0 

-1/2 

0 

2 

N2  4  N2 

T” 

2N  4  N2 

2.  8  x  1010 

-1/2 

0 

2 

N2  +  N 

r- 

2N  4  N 

2.  4  x  1015 

-3/2 

0 

2 

N2  +  NO 

2N  4  NO 

1. 1  x  1010 

-1/2 

0 

2 

NO  4  02 

N  4  O  4  O2 

1. 0  x  10l4 

-3/2 

0 

2 

NO  4  O 

N  4  O  4  O 

1.  0  x  1014 

-3/2 

0 

2 

NO  4  N2 

T— 

N  4  O  4  N2 

1. 0  x  1014 

-3/2 

0 

2 

NO  4  N 

N  4  O  4  N 

1.0  x  1014 

-3/2 

.p 

1 

2 

NO  4  NO 

N  4  O  4  NO 

2.  0  x  1015 

-3/2 

2 

O  4  N2 

NO  4  N 

1.  6  x  1010 

0 

B 

NO  4  O 

N  4  02 

1.  3  x  107 

1 

B 

N2  +  O2 

2NO 

2.  4  x  1020 

-5/2 

mm 

B 

Equilibrium  Reverse  Reaction  Rate  Constant 

^f^EQ  =  |^« 

rC2_/cAl  /  me«erJ  \n 

P  \T  /  J  \  kg-mole  / 
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APPENDIX  III 

TEMPERATURE  MEASUREMENTS  BY  THE  ELECTRON-BEAM  TECHNIQUE 


If  a  low  density  gas  is  bombarded  by  a  high  energy  (10-  to  30-kv) 
beam  of  electrons,  the  resulting  emission  can  be  analyzed  to  infer  the 
gas  temperature.  To  date,  the  only  gases  studied  by  this  technique 
have  been  N2  or  N 2 -rich  mixtures,  using  the  analysis  procedure 
developed  by  Muntz  (Ref.  15).  The  nitrogenoion  (N24)  first  negative 
emission  system,  roughly  centered  at  3900  A,  is  very  intense  in  the 
emission  of  air  excited  by  an  electron  beam.  Muntz  proposed  a  simple 
excitation /emission  model  for  this  emission  system  in  which  the  excited 
N2+  ions  forming  the  upper  state  of  the  emission  were  produced  by 
direct  excitation  of  ground  state  neutral  N2  molecules  in  collision  with 
the  high-speed  electrons  of  the  beam.  The  excited  ions  were  assumed 
to  undergo  spontaneous  emission  to  the  ionic  ground  state  to  produce 
the  first  negative  system,  and  both  excitation  and  emission  were 
assumed  to  obey  optical  selection  rules.  The  three  molecular  states 
involved  have  the  following  designations: 

N2  ground  state  (N2X*E),  designated  by  a  subscript  g 

N24  excited  state  (N24B2£),  designated  by  a  single  superscript 

N24  ground  state  (N24X2£),  designated  by  a  double  superscript 

ROTATIONAL  TEMPERATURE 

In  most  practical  cases,  the  translational  temperature  of  molecular 
motion  is  in  equilibrium  with  the  temperature  characterizing  the  distribu¬ 
tion  of  rotational  energy  states.  The  rotational  temperature  can  be  in¬ 
ferred  from  the  electron-beam  emission  by  analysis  of  the  distribution 
of  intensity  of  individual  rotational  lines  within  a  single  vibrational  band. 
Based  on  the  model  described  above,  Muntz  derived  an  expression  for 
the  intensity  of  a  single  rotational  line: 

Ik  K",  v',  v"  =  constant  x  N0  x  f(v',  v",  Tr,  Tv)  x  F*k K  ",  v',  v" 

x  G(K  0rotf  Tr)  x  K '  exp  |j  K'(K'T  1)]  (m-1) 

where  Ij^  k'',  v',  v''  *s  the  intensity  of  emission  resulting  from  a  rota¬ 
tional  state  transition  from  K'  to  K"  levels,  within  a  vibrational  state 
transition  from  v'  to  v"  level.  Since  the  analysis  is  confined  to  a 
specific  vibrational  band  (usually  0,  0  band),  the  v'  and  v''  subscripts 
may  be  dropped.  Since  optical  transition  rules  are  adopted,  AK  =  ±1  in 
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the  transition,  and  Eq.  Ill- 1  has  already  been  specialized  to  the  case 
where  K''  =  K'  -  1,  the  so-called  R-Branch  of  the  rotational  transition. 
The  product  of  K'  and  the  exponential  function  in  Eq.  ( III- 1 )  is  almost 
identical  to  the  rotational  state  distribution  of  the  unexcited  N2  mole¬ 
cules,  but  in  the  formation  of  the  intensity  distribution  it  is  modified 
by  the  function. 


G(K',  0rot,  Tr)  = 


K'-  1 
2K'+  1 


exp 


—rot  ■  2(K '+  1) 
Tr 


K' 


2K  '+  1 


exp 


— — —  .  2K' 
_  Tn 


(m-2) 


The  determination  of  Tr  from  measured  Ir'  in  a  band  involves  the  re¬ 
arrangement  of  Eq.  ( III  - 1 )  and  the  taking  of  a  natural  logarithm: 


l°ge 


iK'  " 

-  -(  ^  ) 

.fK'GKl 

\  Tr  / 

K'(K'  +  1)  +  constant 


(IH-3) 


Clearly,  there  is  a  linear  relation  between  the  logarithm  on  the  left 
and  the  quantity  K'(K'  +  1)  so  that  a  plot  of  the  logarithm  versus 
K'(K'  4-  1)  will  be  a  straight  line  (see  below)  with  a  slope  equal  to 
-0rot/TR»  where  0rot  for  ground  state  N2  at  Vg  =  0  is  2.  89°K.  The 
"constant"  on  the  right  side  of  Eq.  (III-3)  is  a  function  of  many  vari¬ 
ables:  electron -beam  characteristics,  N0,  the  specific  vibrational 
transition,  Tv  and  Tr.  This  constant  determines  the  intercept  on  the 
ordinate  but  has  no  effect  on  the  slope  of  the  line.  The  function,  G,  is 
weakly  dependent  on  Tr,  necessitating  an  iteration  to  obtain  a  con¬ 
sistent  value  of  Tr.  There  is,  for  this  particular  emission  system, 
a  2:1  intensity  ratio  for  alternating  lines,  resulting  from  nuclear 
statistics.  This  phenomenon  is  usually  compensated  for  by  multiply¬ 
ing  the  half-strength  lines  by  a  factor  of  2.  0. 


i°ge 
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This  procedure,  while  simple  and  straightforward  in  theory,  has 
proved  to  be  somewhat  less  than  simple  to  apply  in  practice.  Muntz 
applied  the  technique  to  measure  static  temperature  in  both  arc -heated 
and  cold  low  density  flows,  and  it  was  rapidly  adopted  by  other  labora¬ 
tories  (Refs.  9,  10,  12,  16,  and  17).  Several  anomalies  were  dis¬ 
covered  when  the  technique  was  extended  to  temperature  regimes  not 
covered  in  Muntz's  original  experiment.  First  of  all,  at  temperatures 
below  150°K  the  logarithmic  intensity  plots  were  found  to  exhibit  a  pro¬ 
nounced  upward  curvature  which  increased  as  absolute  temperature  was 
decreased.  Second,  even  when  straight  lines  were  drawn  tangent  to 
these  curves,  the  indicated  Tr  was  invariably  in  excess  of  the  gas  tem¬ 
perature  obtained  by  gas  dynamics  calculations,  and  the  excess  also 
increased  as  temperature  decreased.  Since  these  effects  were  first 
observed  in  low  temperature  free  jets  and  nozzle  flows  (Ref.  16),  the 
possibility  of  rotational  nonequilibrium  in  the  flow  expansions  was  sug¬ 
gested.  This  possibility  was  eliminated,  however,  by  an  extensive 
experiment  conducted  in  a  static  test  chamber  at  temperatures  down 
to  80°K  (Ref.  20).  The  same  results  were  found  as  for  a  flowing  gas, 
except  that  an  additional  dependence  of  the  divergence  of  Tr  and  gas 
temperature  on  static  density  was  found.  A  number  of  possible  explana¬ 
tions  of  these  effects  has  been  investigated,  including  a  complete  re¬ 
formulation  of  the  excitation  model  (Ref.  38),  but  their  source  still 
remains  uncertain. 

At  very  low  temperatures  (10  to  20°K)  the  divergence  of  Tr  and  gas 
temperature  is  as  large  as  40  percent,  but  at  the  temperatures  normally 
encountered  in  the  18-in.  low  density  wind  tunnel  (100°K),  it  is  only 
approximately  8  percent.  The  measured  temperatures  are  easily  cor¬ 
rected  by  use  of  an  empirical  correction  curve,  as  suggested  in  Ref.  16. 
The  density  dependence  of  the  divergence  is  rather  small  (5  percent  in 
an  order  of  magnitude),  so  that  it  is  usually  possible  to  use  a  single 
correction  curve  over  the  entire  range  of  operation  of  a  given  facility. 

In  the  present  study,  the  correction  curve  of  Ref.  16  was  used  directly, 
except  for  minor  changes  made  to  it  at  temperatures  around  80°K, 
based  on  static  calibrations  in  our  own  facilities.  The  wind  tunnel  test 
results  have  shown  that  the  correction  is  indeed  required  and  that  the 
empirical  correction  of  Ref.  16  appears  to  be  of  the  proper  magnitude. 


VIBRATIONAL  TEMPERATURE 

The  temperature  characterizing  the  distribution  of  vibrational  states 
of  the  unexcited  N2  molecules  is  obtained  by  reference  to  the  total  in¬ 
tensity  (sum  of  all  rotational  lines)  of  a  specific  vibrational  band  in  the 
N2+  emission.  Using  Muntz's  excitation /emission  model,  it  can  be 
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shown  that  the  intensity  of  a  complete  vibrational  band  representing  a 
transition  between  vibrational  levels  specified  by  the  quantum  num¬ 
bers  v'  and  v''  is : 


=  constant  x  N0  x 


v4vV'p(v',  v\S=fxp  [-  4^-]  Pb'.Vg) 


Qv(Tv) 
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The  ratio  of  intensity  of  two  different  bands  is  thus. 
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For  any  two  specific  bands,  I1/I2  =  f(Tv),  a  function  which  has  the 
general  shape  shown  below.  For  the  present  test,  this  curve  was  com¬ 
puted  for  the  (0,  1) /<  1,  2)  band  ratio  and  used  to  determine  the  Tv  given 
in  Table  I. 


Unfortunately,  the  transition  probability  factors  for  the  vibrational 
transitions  are  not  known  with  the  same  degree  of  accuracy  as  are  the 
analogous  factors  for  the  rotational  transitions,  and  there  is  consider¬ 
able  difference  of  opinion  over  the  accuracy  of  Tv  obtained  by  this  pro¬ 
cedure.  Muntz  originally  quoted  an  accuracy  of  ±5  percent,  but  later 
investigators  have  been  less  generous  in  their  estimates,  with  values 
ranging  from  ±10  to  ±20  percent  (Refs.  9  and  10).  Further,  very  little 
has  been  done  to  experimentally  verify  the  measurements  under  the 
electron-beam  conditions.  A  static  calibration  experiment  is  reported 
in  Ref.  39  in  which  three  points  between  300  and  1100°K  were  obtained. 
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showing  a  divergence  of  from  1.  5  to  16  percent  from  the  theoretical 
curve.  In  Ref.  10,  Tv  obtained  in  Wind  tunnel  flows  using  three  dif¬ 
ferent  band  ratios  were  found  to  have  a  scatter  of  ±7  percent  between 
the  three  ratios  used.  In  the  present  study,  a  single  comparison  of  Tv 
from  two  different  band  ratios  showed  a  similar  7-percent  difference. 

A  recent  study  of  the  electron-beam  temperature  measurement  tech¬ 
nique,  giving  detailed  derivations  and  results  of  validation  experiments, 
is  reported  in  Ref.  40. 
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APPENDIX  IV 

DETERMINATION  OF  DENSITY  BY  ELECTRON  BEAM 


There  are  two  types  of  electron-beam  techniques  which  have  suffi¬ 
cient  accuracy  and  resolution  to  be  used  for  wind  tunnel  diagnostics. 

The  first  of  these  utilizes  the  emission  of  the  target  gas  as  a  measure 
of  gas  density.  In  Eqs.  <  III  - 1 )  and  (III-4),  it  can  be  seen  that  both  indi¬ 
vidual  rotational  lines  and  complete  vibrational  bands  have  intensities 
directly  proportional  to  No,  the  number  density  of  neutral  N2  molecules. 
The  various  applications  of  this  technique  differ  in  the  emission  com¬ 
ponents  used.  Muntz  used  the  intensity  of  a  single  rotational  line  in  the 
(0,  0)  band  (Ref.  15).  In  Ref.  9  the  total  (0,  0)  band  strength  was  used, 
and  in  Ref.  10  the  (0,  2).  band  was  used.  In  Ref.  17,  a  number  of  bands 

O 

were  used  simultaneously  by  means  of  a  100  A  half-width  interference 
filter  centered  at  the  (0,  0)  band.  This  technique  is  very  easily  cali¬ 
brated  by  measurements  in  a  static  gas  at  precise  values  of  pressure 
and  temperature.  If,  however,  the  rotational  and  vibrational  tempera¬ 
tures  prevailing  during  the  flow  measurement  are  not  the  same  as  in 
the  static  calibration,  correction  factors  are  required  which  depend  on 
these  temperatures.  This  technique  is  also  affected  by  any  chemical 
changes  in  the  gas  which  change  the  relative  N2  concentration. 


The  second  type  of  density  measurement,  and  the  one  used  in  the 
present  study,  is  based  on  measurement  of  electrons  elastically 
scattered  at  large  angles  out  of  the  beam  volume.  A  comprehensive 
description  of  the  development  of  this  technique  is  given  in  Ref.  19. 
The  differential  scattering  cross  section  which  an  atomic  nucleus  pre¬ 
sents  to  a  high-speed  electron,  which  penetrates  the  atom  electron 
cloud  and  approaches  the  nucleus,  is  (Ref.  19) 


d  a 


(0  sin  -|-) 


dQ 


(IV-1) 


where  the  deflection  is  through  an  angle  6  and  into  a  solid  angle  of  dfl. 
If  an  electron  beam  is  directed  into  an  extended  region  of  target  atoms 
of  particle  density  n,  there  will  result  a  field  of  scattered  high-speed 
electrons  surrounding  the  beam,  with 
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Z2  (1  -  ff2) 
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where  X  is  the  rate  at  which  electrons  are  scattered  from  a  small  seg¬ 
ment  of  the  beam  Ax  into  a  solid  angle  AQ  at  an  angle  6  with  respect  to 
the  incident  beam.  Rearranging  this  expression  into  convenient  groups 
of  variables. 
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The  beam  characteristics  are  determined  by  current  and  voltage,  and 
the  detector /collimator  characteristics  are  mainly  a  matter  of  geo¬ 
metrical  design.  The  target  gas  characteristic  enters  only  as  the 
atomic  number  of  the  nucleus  and  the  particle  density  of  the  nuclei.  In 
this  model  of  the  scattering  process,  based  on  deep  penetration  of  the 
atom  electron  cloud,  the  distribution  of  nuclei  through  various  chemical 
combinations  does  not  affect  the  scattering  process  to  first  order.  The 
mass  density  of  a  mixture  can  be  expressed  as 


N, 

where  the  summation  depends  only  on  the  total  number  of  different 
types  of  atoms  present  in  the  mixture.  Thus,  a  calibration  of  X  versus  p 
made  in  low  temperature  air  will  still  be  valid  at  high  temperature  even 
though  appreciable  changes  in  detailed  composition  may  have  occurred. 
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